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ABSTRACT 

Context. We present detection and analysis of faint X-ray sources in the Chandra deep field south (CDFS) using the 4 
Msec Chandra observation. 

Aims. We put new constraints on the z=3-7 AGN luminosity function, its cosmological evolution and high-redshift 
black hole and AGN demography. 

Methods. We use a new detection algorithm, using the entire 3-dimensional data-cube (position and energy), and use a 
targeted search at the position of high-z galaxies in the GOODS-South survey. 

Results. This optimized technique results in the identification of 54 z>3 AGNs, 29 of which are new detections. Applying 
stringent completeness criteria, we derive AGN luminosity functions in the redshift bins 3-4, 4-5 and >5.8 and for 
42.75<log L(2-10 keV)<44.5. We join this data with the luminous AGN luminosity functions from optical surveys 
and find that the evolution of the high-z, wide luminosity range luminosity function can be best modeled by pure 
luminosity evolution with L* decreasing from 6.6 x 10*^ ergs/s at z=3 to L* = 2 x 10^* ergs/s at z=6. We compare 
the high-z luminosity function with the prediction of theoretical models using galaxy interactions as AGN triggering 
mechanism. We find that these models are broadly able to reproduce the high-z AGN luminosity functions. A better 
agreement is found assuming a minimum Dark Matter halo mass for black hole formation and growth. We compare our 
AGN luminosity functions with galaxy mass functions to derive high-z AGN duty cycle using observed Eddington ratio 
distributions to derive black hole masses. We find that the duty cycle increases with galaxy stellar mass and redshift by 
a factor 10-30 from z=0.25 to z=4-5. We also report on the detection of a large fraction of highly obscured, Compton 
thick AGN at z> 3 (18l]^g%). Their optical counterparts do not show any reddening and we thus conclude that the 
size of the X-ray absorber is likely smaller than the dust sublimation radius. We finally report the discovery of a highly 
star-forming galaxy at z=3.47, arguing that its X-ray luminosity is likely dominated by stellar sources. If confirmed, 
this would be one of the farthest objects in which stellar sources are detected in X-rays. 

Key words. Galaxies: active - Galaxies: evolution - Galaxies: quasars - general 

1. Introduction IBoutsia et al. 20TT1) . However, there are other funda- 
mental issues that can be tackled by studying high-z 

The study of high redshift AGN holds the key to under- ^GN: a) scenarios for the formation of the black hole 

stand early structure formation and probe the Universe (^H) seeds, which will eventually grow up to form the 

during its infancy High-z AGN have been extensively super-massive black holes (SMBHs) seen m most galaxy 

used to investigate key issues such as the evolution bulges. Two main scenarios have been currently proposed, 

of the correlations between t he black hole mass and the monoUthic collapse of 1 - 10^ M0 gas clouds to BH 

the galaxy properties (see e.g ILamastra et al. 20101 and dLodato fc Natara.jan 2006[ |Lodato fc Natarajan 2007[ 



references therein), the AGN contribution to the re- |Begelman 201 0[ [Voionteri fc Begelman 2010p , or e^ 

ionization, heating of the Inter- Ga lactic Matter and its generation of ~ 100 M© BH prod uced by the supernov a 

effect on structure formation (e.g. [Giallongo et al.T997l explosions of the first Pop HI stars (|Madau fc Rees 200 ip. 

IMalkan et al. 20031 [Hopkins et al. 2007[ IGowie et al. 20101 b) Investigate the physics of accretion at high-z: is 
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BH growth mainly due to relatively few accretion 
episodes, as predicted in hierarchical scenarios (see e.g. 
[Dotti et al. 2010 and references therein), or by the so 
called chaotic accretion (hundreds to thousands of small 
accretion episodes, King et al. 2008[ )? The two scenarios 
predict different BH spin distributions, and thus different 
distributions for the radiative efficiency, c) Being BH the 
structures with the fastest (exponential) growth rate, 
they can be used to constrain both the expansion rate 
of the Universe and the growth rate of the primordial 
perturbation at high-z, i.e. competing cosmological sce- 
narios (jFiore 20101 ILamastra et al. 2011|) . d) The slope of 
the high-z AGN luminosity function and of the SMBH 
mass function strongly depend on the AGN duty cycle, 
and therefore their measure can constrain this critical 
parameter. In turn, the AGN duty cycle holds information 
on the AGN triggering mechanisms. The evaluation of the 
evolution of the AGN duty cycle can thus help in disen- 
tangling among competing scenarios for AGN triggering 
and feeding. In this paper we discuss in more detail this 
last issue. Two main scenarios for AGN triggering and 
feeding have been investig ate d so far: galaxy encounters 
(Barnes fc Hernquist 1996 

and recycled gas 



ICavahere fc Vittorini 20001 

IMenci et al. 20061 IMenci et al. 200"8| . 



produced by normal stellar evolution in the inner bulge 
(ICiotti fc Ostriker 2003 IGiotti et al. 20T0l lUen 2011|) . 
These make different predictions about the AGN time 
scale and duty cycle. In the former models nuclear activa- 
tion follows galaxy interactions. Since these are naturally 
more frequent in the past, when, at the same time, more 
gas is available for nuclear accretion since it has not yet 
been locked in stars, this scenario predicts a strong increase 
of the AGN duty cycle with the redshift. In the recycled 
star gas scenario the AGN timescale is much longer, up 
to Gyr, although with decreasing Eddington ratio. A less 
extreme variation of the AGN duty cycle with the redshift 
is thus expected. One of the goals of this paper is to assess 
what can be already done with present AGN surveys to 
distinguish between these competing scenarios, and which 
kind of surveys should be planned to best assess the main 
AGN triggering and feeding mechanism. 

Large area optical surveys such as the SDSS, 
the CFHQS, and the NOAO DWFS/DLS have al- 
ready been able to discover large samples of z> 4.5 
QSOs (e.g. [Richards et al. 20061 IGlikman et al. 201 ip and 
about 50 QSOs a t z > 5.8 (e.g. [Jiang et al. ^009 
IWillott et al. 2010a|) . The majority of theseTiighTTGI^ 
are broad line, unobscured, high luminosity AGN. They 
are likely the tips of the iceberg of the high-z AGN popula- 
tion. Lower luminosity and/or moderately obscured AGN 
can, in principle, be detected directly in current and future 
X-ray surveys. Dedicated searches for high-z AGN using 
both deep and wide area X-ray surveys and a multi-band 
selection of suitable candidates can increase the number 
of high-z AGN by a factor > 10. In particular, it should 
be possible to find hundreds rare high-z, high luminosity 
QSOs, in both the all sky and deep eROSITA surveys (the 
0.5-2 keV flux limit of the all sky survey being the order of 
10"""^^ erg/cm^/s , while that of the deep survey, covering 
hundreds deg^, should be 2-3 times deeper) with a selec- 
tion function much less biased than optical surveys (recall 
that Tx/Topt w (1 -I- z)~^'^, where tx and Topt are the opti- 
cal depths in the observed- frame X-ray and optical bands, 
assuming no evolution of the dust-to-gas ratio). To con- 



strain the faint end of the high-z AGN luminosity func- 
tion, and therefore the shape of the luminosity function 
and of the SMBH mass function, we need to best exploit 
current and future deep surveys. Unfortunately, so far the 
number of X-ray selected AGN at z> 3 in deep Chandra 
and XMM surveys is only of a few tens, and ~half a dozen 
at z> 4.5 (see e.g. Brandt et al. 2004, Fontanot et al. 20071 
,Brusa et al. 2009 a, Brusa et al. 2009b, Civano et al. 201 ip . 
The difficulty in detecting directly high-z AGN in X-ray 
surveys suggests an alternative approach. Looking at the 
X-ray emission at the position of known sources allows 
one to use a less conservative threshold for source detec- 
tion than in a blind search, and therefore to reach a lower 
flux limit. Furthermore, a highly optimized X-ray analy- 
sis tool is deemed mandatory to fully exploit the richness 
of the data produced by Chandra and XMM. To this pur- 
pose we developed a new X-ray detection and photome- 
try tool (dubbed ephot, see the Appendix). In this paper 
we search for X-ray emission at the position of candidate 
high-z galaxies selected in the HST/WFC3 ERS (Early re- 
lease science) area and in the GOODS area of the Chandra 
Deep Field South (CDFS) by using ephot. We use the new 4 
Msec Chandra dataset, which represent the most sensitive 
X-ray exposure ever achieved. We first study selection ef- 
fects affecting our candidate high-z AGN, and then join our 
CDFS samples to other X-ray and optically selected AGN 
samples at the same redshift to build AGN luminosity func- 
tions at z> 3 over several luminosity decades. By assuming 
appropriate bolometric corrections and appropriate distri- 
butions of A = Liioi I L Edd we convert these luminosity func- 
tions in SMBH mass functions. We then evaluate the stellar 
mass functions of active galaxies by converting the SMBH 
mass to galaxy stellar mass by using the Magorrian rela- 
tionship (|Ferrarese fc Ford 20051 IShankar et al. 2009al and 
references therein). Finally we estimate the AGN duty cy- 
cle as a function of the redshift by dividing the stellar mass 
function of active galaxies by the stellar mass function of 
all galaxies at the same redshift. The AGN duty cycle as a 
function of the redshift is then compared with the expecta- 
tion of models for structure formation using galaxy inter- 
action as AGN triggering mechanism. K Hq — 7Q km s~^ 
Mpc~^, r2j\/=0.3, 17a — 0.7 cosmology is adopted through- 
out. 



2. Candidate high-z galaxy samples 

We built a catalog of z> 3 galaxy candidates in the GOODS 
south and ERS fields including: 

— all galaxies with high quality spectroscopic redshift > 3; 

— all galaxies with 68% limit of photometric redshift > 3; 

— all galaxies with B-V vs V-I or V-H colors suggesting 
3.5<z<4.4; 

— all galaxies with V-I vs I-Z or I-H colors suggesting 
4.4<z<6.0; 

— all galaxies with I-Z vs Z-Y or Z-Y vs Y-J colors sug- 
gesting z>6; 

— all galaxies with Y-J vs J-H colors suggesting z>8. 



2.1. GOODS field 



We use the GOOD S-MUSIC catalog (IGrazian et al. 20061 
ISantini et al. 2009|) . based on the mosaics reduced as de- 
scribed in Koekemoer et al. (2011), see also Windhorst et al. 
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(2011), which includes galaxies selected in the HST/ACS z 
band and Spitzer/IRAC 4.5//m band. 190 GOODS-MUSIC 
galaxies have a high quality optical or near-infrared spec- 
trum confirming that they are z> 3 galaxies, 26 have a 
spectroscopic confirmation at z> 5. 4284 galaxies without 
high quality spectroscopic redshift •L^-pec have photometric 
redshift i^hot > 3, considering the 68% upper limit. To 
this sample, we added 171 galaxies without high quality 
Zspec and Zphot < 3 that lie in the 3.5<z<4.4 region of the 
ACS B-V vs V-I diagram and 143 galaxies that lie in the 
4.4<z<6.0 region of the ACS V-I vs I-Z diagram. The total 
number of high-z galaxy candidates considered is thus 4788. 
About 10% of these galaxies lie within a few arcsec from 
an X-ray source previously identified with a low-z AGN or 
galaxy or galaxy cluster. These objects were conservatively 
excluded from the list, leaving us with about 4300 high-z 
galaxy candidates. 

2.2. ERS field 

We use the GOO DS-ERS catalog (|Grazian et al. 20101 
ISantini et al. 20ll|l . which includes 2291 galaxies selected 
in the HST/WFC3 H band and with photometric redshift 
>3. 42 of these objects have a confirmed spectroscopic red- 
shift (only 2 with z> 5). Again, about 10% of these galaxies 
lie close to bright X-ray sources, thus limiting the ERS can- 
didate high-z galaxy sample to about 2000 galaxies. 

There is of course an overlap between the GOODS-ERS 
and the GOODS-MUSIC catalogs (the GOODS field covers 
nearly completely the ERS field). Of the 2291 GOODS- 
ERS, z> 3 galaxies 570 are present in the GOODS-MUSIC 
catalog and in 420 cases the GOODS-MUSIC photometric 
redshift is > 3 Therefore, these objects are present also in 
the GOODS-MUSIC high-z galaxy candidate sample. 

2.3. Additional samples 

We also considered 330 z> 6 galaxy candi- 
dates from the lists of IBouwens et al. 20"06l 

ICastellano et al. 20101 IMcLure et al. 20101 

IMcLure et al. 20Tn IBouwens et al. 20101 which are 
not included in the above two lists. 

3. Multiwavelength properties of High-z AGN and 
galaxies in the CDFS. 

To extract statistically quantitative information from a 
given sample of "detected" sources in an X-ray image, 
we need to assess its reliability (i.e. number of spurious 
detections) and completeness. To both purposes we 
run a series of extensive detection runs on simulated 
data, see the Appendix for details. We choose a detec- 
tion threshold corresponding to 1 spurious detection 
in 5000 candidates, i.e. about 1 spurious detection in 
the GOODS-MUSIC sample and <1 spurious detection 
in the GOODS-ERS sample. We detect 17 sources at 
z> 3 from the GOODS-ERS catalog and 41 from the 
GOODS-MUSIC catalog within this threshold. As an 
example. Fig. [T] show the B, z and H band images of 
three Chandra-GOODS-ERS sources. Table 1 and 2 give 
position, redshift. X-ray band that optimizes the detection, 
0.5-2 keV flux, z band and H band magnitudes and 4.5/im 
fluxes for the GOODS-ERS and GOODS-MUSIC X-ray 
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Fig. 1. HST ACS B band (left panel), z band (central 
panel) and HST WFC3 H band (right panel) images of 
three Chandra-GOODS-ERS sources with overlaid X-ray 
contours in the 0.7-2.7 keV band (E2199), 0.7-1.9 keV band 
(E2551) and 0.4-1.5 keV band (E4956) 



detected sources, along with a few sources just below 
detection threshold (4 from the ERS-GOODS sample, 
2 from the GOODS-MUSIC). We do not detect any 
of the candidate z ^ 6 galaxies in IBouwens et al. 20"06l 

ICastellano e t al. 20101 IMcLure et al. 20101 

McLure et al . 20111 IBouwens et al. 20101 One of these 
galaxies has an X-ray counterpart just below threshold 
(also reported in Table 2). At the position of this source 
there are 11 counts in the band 0.4-0.8 keV 3 of which 
are background counts) in a 2 arcsec radius region. If the 
0.4-0.8 keV emission is due to iron Ka emission at 6.4 keV 
the redshift would be ^8, which is inconsistent with the 
detection in the z band. Because of these uncertainties we 
do not discuss further this source. 

Five Chandra-GOODS-ERS sources are present in the 
2 Msec direct detection list of iLuo et al. 20081 (Luo08 here- 
after, search within a circle of 1 arcsec radius). Three of 
them have Zphot in ILuo et al. 20101 (LuolO hereafter) con- 
sistent with those in the GOODS-ERS catalog, while in 
one case the LuolO Zphot is just below 3. One Chandra- 
GOODS-ERS source is confused/blended in the Luo08 cat- 
alog. Eight Chandra-GOODS-ER S sources are pre sent in 
the 4 Msec direct detection list of IXue et al. 20lT] fXuell 
hereafter). Two of the three Xuell sources not in the Luo08 
catalog have a photometric redshift, and this is > 3 in both 
cases. 16 Chandra-GOODS-MUSIC sources are present in 
the LuoOS list, 13 with photometric or spectroscopic red- 
shift consistent with those in the GOODS-MUSIC catalog. 
Three additional sources are present in the Xuell catalog, 
two of them have a photometric redshift, and this is con- 
sistent with that in the GOODS-MUSIC catalog. Source 
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M208 was not detected by Luo08, while it is now detected 
by Xuell. It is extensively discussed in lGilli et al. 20lT] 

The accuracy of the photometric redshifts is greatly 
improved thanks to the near infrared bands covered by 
WFC3 (see Table 1 and 2). The WFC3 coverage also in- 
creases the high-z candidate galaxy X-ray detection rate. 
The ERS area is less than one third of the fuU GOODS 
area, and therefore, based on the Chandra-GOODS-ERS 
detections, one would expect detections in GOODS- 
MUSIC, while only 41 have been actually found. Only half 
of the Chandra-GOODS-ERS detections have a counter- 
part in the GOODS-MUSIC catalog. 31% and 40% of the 
Chandra-GOODS-MUSIC detections are also present in the 
LuoOS 2 Msec and Xuell catalogs. 28% and 47% of the 
Chandra-GOODS-ERS detections are present in the same 
catalogs. We study in more detail the different Chandra- 
GOODS-MUSIC and Chandra-GOODS-ERS selection ef- 
fects in section 3.3. 

3.1. X-ray and optical spectroscopy 

We discuss here X-ray and optical spectroscopy of our X- 
ray sources. Unfortunately only one of the 17 sources in 
the Chandra-GOODS-ERS sample has optical spectroscopy 
available. On the other hand, 11 of the 41 sources of 
the Chandra-GOODS-MUSIC sample have optical spec- 
troscopy available. To exploit this rather large dataset we 
complement the full Chandra-GOODS-ERS sample, with 
the Chandra-GOODS-MUSIC sources with optical spec- 
troscopjQ. 

X-ray spectra have been extracted from the merged 
event file using source extraction regions of 2-6 arcsec 
radius, depending on the source off-axis (the Chandra 
Point Spread Function degrades quickly with the off- 
axis angle). Background spectra have been extracted from 
nearby source-free regions. We also used the total back- 
ground spectra at three off-axis angles (see the Appendix 
for details). Background spectra are then normalized to 
the source extraction region area. We report in Table 
3 the results of fitting a simple power law model to 
the data. Several Chandra-GOODS-ERS and Chandra- 
GOODS-MUSIC sources have flat X-ray spectra, strongly 
suggesting the presence of large absorbing column densi- 
ties along those lines of sight. The quality of the spectra is 
poor, however in several cases the statistic is good enough 
to perform a proper fit with a model including an absorbed 
power law or a reflection spectrum. In the following we 
briefly outline the main results from this analysis. 

— E537 shows an extremely hard spectrum. A simple 
power law fit produces a negative slope of F ^ -0.15. 
An absorbed power-law model with F fixed to 1.8 yields 
a column density of N// = 1.6 x lO^"' cm~^. A better de- 
scription of the spectral shape of the low energy spec- 
trum is obtained by adding an unabsorbed power law 
(see Fig. [2]) . Such component typically provides an ex- 
cellent fit of the X-ray spectrum of many well-studied 
obscured AGNs and its origin is explained in terms 
of emission due to reprocessing of the nuclear contin- 

^ We do not discuss here the rest of the Chandra-GOODS- 
MUSIC sample because it is affected by incompleteness, see next 
section, and because the quality of the photometric redshift is 
worse than that of the GOODS-ERS sample 



uum by surrounding material and/or a portion of pri- 
mary radiation leaking through the absorber (see e.g. 
[Turner et aTT 997, Piconcclh ct al. 2007). The photon 
index of the soft X-ray component is fixed to F=1.8, 
while the resulting ratio between the normalization of 
the unabsorbed and absorbed power law is ~0.02. An al- 
ternative model consisting of a Compton refiection emis- 
sion component (pexrav model in XSPEC) obscured by 
a column density of — 5±2 x 10^^ cm^-^ is statisti- 
cally equivalent. 

— The X-ray spectrum of M4835 reveals additional com- 
plexity when fitted with a simple absorbed (Nh = 9±2 
X 10^^ cm"2) power-law model (C-stat/dof=425/442). 
The best-fit model to the Chandra data includes an ad- 
ditional component due to X-ray reflection from cold 
circumnuclear matter. In this case the primary contin- 
uum is obscured by a column density of Nh ~ l-5to 8 ^ 
lO'^'* cm~^ and the normalization of the nuclear compo- 
nent is about 5 times larger than that of the reflection 
component. Finally, the limit equivalent width (EW) of 
the Fe Ka emission line measured with respect to the re- 
flection continuum is <1.1 keV, i.e. still consistent with 
the expectation of a reflection model. Fig. [5] shows the 
Chandra spectrum of M4835 fltted by this reflection -t- 
transmission model. 

— M5390 is G202 the well-known, hi ghly obscured QSO 
discovered by INorman et al. 20021 at z = 3.7. The 
Chandra spectrum is consistent with the XMM one pre- 
sented by IComastri et al. 20111 It is well fitted by a 
power law model (with F = 1.8 fix) absorbed by a col- 
umn density ofN// '-^ 1.2±0.15xl0^'' cm~^ and gaussian 
line (E=6.6 keV, likely a blend of the 6.4 keV neutral 
Fe and 6.7 keV Fe XXV Ka hues, EW=0.8±0.4 keV), 
plus an energetically unimportant soft component. The 
quality of the fit decreases if the continuum is modeled 
with a pure reflection component. 

— M208 is the source presented bv lGilli et al. 20111 and we 
confirm the results reported in that paper, finding the 
spectrum consistent with being obscured by a Compton 
thick absorber. Very similar results are obtained for 
E1577, which is best fitted by a heavily absorbed (i.e. 
N_f/ > 10^** cm~^) power law model (see Fig. We 
also find the presence of obscuration with column den- 
sity values consistent with 10^^ cm~^ in the spectrum 
of E8479. 

— Large N// values exceeding 10^^ cm~^ (but still in 
the Compton-thin regime) are reported from the spec- 
tral analysis of E1617, M8273, M4302 and M3320. 
Remarkably, in the spectrum of the latter source at 
z=3.471, we find clear evidence of an emission line with 
an EW — 160±100 eV at a best-fit rest-frame energy of 
6.97 keV and, therefore, associated with highly ionized 
(i.e. Fe XXVI) iron. 

In conclusion, at least 3 of the 17 Chandra-GOODS- 
ERS sources are likely to be Compton thick AGN, which 
is 18l|^o^ AGN in the sample are Compton thick, 

as weU as 4 of the 11 Chandra-GOODS-MUSIC AGN 
with spectroscopic redshift (one of these sources is in 
common with the Chandra-GOODS-ERS sample). Other 
3 Chandra-GOODS-MUSIC and 1 Chandra-GOODS-ERS 
AGN are highly obscured. The luminosity of the Compton 
thick ERS and GOODS-MUSIC AGN' is in the range 
log(L2-10keV)=43. 5-44.5, i.e. between bright Seyfert 2 
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Table 1. Chandra detected GOODS-ERS z> 3 galaxies 
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HA 


Dec 
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zm 


TT 
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E537 


53.082535 


-27.755245 


4.29 (2.9-4.5) 


1.3-6.3 


13.6 


1.2E-16 


25.48 


24.21 


L184, X275, M70139'' 


E737 


53.050343 


-27.752138 


4.14 (4.0-4.4) 


1.0-2.0 


3.4 


4.5E-17 


25.65 


24.94 


X177, M70131'' 


E1312 


53.016602 


-27.744848 


3.48 (3.4-3.6) 


0.9-2.1 


5.6 


l.OE-16 


25.30 


24.39 


L68, XlOO, M13365'', MY5063 


E1516 


53.119888 


-27.743042 


4.29 (4.1-6.9) 


0.5-1.9 


3.4 


3.2E-17 


27.99 


25.57 


X392, M70099 


E1577 


53.020576 


-27.742151 


3.462 spec. 


1.9-6.3 


3.7 


9.5E-17 


23.62 


23.35 


M13549,MY51260 


E1617 


53.010651 


-27.741613 


3.38 (3.3-3.4) 


0.9-3.8 


7.7 


1.9E-16 


26.33 


24.58 


L61^X90 


E2199 


53.054756 


-27.736839 


4.71 (4.5-5.2) 


0.7-2.7 


2.9 


3.5E-17 


27.15 


27.08 




E2309 


53.117855 


-27.734295 


3.61 (3.3-3.8) 


0.3-4.0 


5.3 


4.7E-17 


23.81 


22.93 


L252, X386, M14078 


E2498 


53.006237 


-27.734076 


6.62 (5.9-7.3) 


1.0-5.3 


5.0 


5.1E-17 


-28.89 


26.61 


L57^, X85 


E2551 


53.158340 


-27.733500 


5.19 (5.1-5.5) 


0.7-1.9 


4.3 


5.0E-17 


-27.59 


25.00 


G557,L323,X521,M70091 


E2658 


53.055538 


-27.732950 


3.74 (3.5-4.1) 


0.4-6.0 


3.8 


4.2E-17 


26.84 


26.85 






CO nQ7vnQ 
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97 70077^^ 


O.Oi yo.yj '^.^j 


n 9 n 


o.o 


7 OR 1 7 




9fi Q4 


T 1 nn^ 


E6257 


53.168419 


-27.719418 


3.22 (3.2-3.5) 


0.8-1.9 


3.1 


4.1E-17 


24.64 


24.56 


X554, M15144/M15156'', MY 


E7911 


53.194885 


-27.699949 


4.86 (4.7-5.4) 


0.8-1.8 


3.2 


1.3E-16 


99.00 


25.68 




E8479 


53.167572 


-27.694883 


3.67 (3.3-3.9) 


1.5-3.2 


3.2 


7.0E-17 


26.56 


26.27 


MY58818 


E10247 


53.131218 


-27.685293 


3.95 (3.6-4.0) 


0.5-4.5 


3.1 


3.2E-17 


26.55 


26.59 




E1611'' 


53.079372 


-27.741625 


6.91 (6.0-7.4) 


0.8-1.7 


2.6 


2.5E-17 


26.15 


23.75 


M70107 


E7840'' 


53.173592 


-27.661259 


3.68 (3.3-4.0) 


0.3-6.3 


2.9 


4.9E-16 


25.04 


25.44 


MY64565 


E8439'' 


53.183681 


-27.694748 


4.91 (4.8-5.2) 


0.8-2.5 


2.9 


3.2E-17 


99.00 


25.26 




E9850'' 


53.102806 


-27.679691 


3.93 (3.8-4.1) 


0.9-4.7 


2.9 


2.1E-16 


26.36 


26.46 





''L=Luo08, M=GOODS-MUSIC, MY=GOODS-MUSYC Cardamone et al. 2010;X=Xuell; ''GOODS-MUSIC z=2.28; '^GOODS- 
MUSIC z=2.33; ''no redshift in GOODS-MUSIC catalog; "^LuolO z=2.123; ■'^LuolO z=2.808; sLuoOS confused with optically bright 
galaxy; ''M15144 z=0.15 M15156 z=3.07; *secondary solution at z~2;'' just below detection threshold. 




galaxies and type 2 QSOs. The 2-10 keV flux of the 
Chandra-GOObS-ERS sources is between 0.3 an d 3 x IQ-^s 
erg/cm /s . At this fluxes the IGilli et al. 20071 model for 
the Cosmic X-ray background (CXB) predicts a fraction 
of Compton thick AGN of ~ 20%. This is already similar 
to the fraction of Compton thick AGN we find at z> 3, 
thus suggesting that Compton thick AGN are probably 
more common at high-z than previously predicted (also see 
IGilh et al. 20TT]) . 

It is instructive to compare the result of the X-ray spec- 
troscopy to those of the optical spectroscopy, photometry 
and morphology in Table 3. The two CT AGN for which 
CIV is redshifted in the band covered by optical spec- 
troscopy show a narrow CIV, in addition to narrow Lya. 
One of these (E1577/M13549) also shows broad absorp- 
tion blueward CIV. The highest redshift CT AGN (M208) 
shows a narrow Lya and a relatively broad (2000km/s) NV 
emission (also see IGilli et al. 20TT|) . Two of the three highly 
obscured AGN (M3320 and M8273) also show narrow Lya 



emission and, when redshifted in the band covered by the 
spectroscopy, narrow CIV. In the spectrum of the third 
highly obscured source (M4302) strong and broad absorp- 
tion lines are present. 

Four of the ten highly obscured and Compton thick 
AGN have a point like morphology in the z band 
(M13549/E1577, M5390, M4302, M208) or H band (E1577) 
based on CLASS_STAR parameter of SEXTRACTOR, 
FWHM measurements, and eye-ball inspection of z and 
H band images. This suggests that the active nucleus con- 
tributes significantly to the UV and optical rest-frame emis- 
sion, as further confirmed in M4302 and M13549/E1577 
by the detection of likely nuclear broad absorption lines in 
their UV rest frame spectra. 

We estimated the rest-frame UV (0.16/im) extinction 
of the highly obscured and Compton thick AGN by fitting 
their observed Spectral Energy Distribution with galaxy 
templates and the lCalzetti et al. 20001 extinction law. Table 
3 gives the la upper limit to A(0.16/im) along with the 
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Table 2. Chandra detected GOODS-MUSIC z> 3 galaxies 



ID 


RA 


Dec 


z 


^min — max 


SNR 


F(0.5-2keV) 

erg/cm /s 


zm 


H 


F(4.5^m) 


other ID" 




J2000 


J2000 




keV 








MJy 




1^208" 


53.122036 


-27.938740 


4.762 spec. 


1.1-4.0 


3.6 


l.OE-16 


24.99 


23.42 


3.0 


G618,X403,A 


M387 


53.088505 


-27.933449 


3.130 (2.4-3.5) 


0.5-1.3 


2.7 


6.7E-17 


25.02 


99.00 


1.1 




M603 


53.109451 


-27.928093 


3.360 (2.2-3.9) 


0.8-1.9 


2.9 


3.7E-17 


25.82 


25.09 


0.5 


MY21391 


M2690 


53.163532 


-27.890448 


3.620 (1.5-4.6) 


0.7-6.0 


3.0 


l.lE-16 


25.67 


24.96 


9.1 


MY27345 


M3268 


53.083511 


-27.881855 


3.300 (2.6-3.8) 


1.7-6.3 


3.2 


l.OE-16 


26.46 


26.47 


0.2 


MY28191 


M3320 


53.209370 


-27.881086 


3.470 spec 


0.7-5.0 


41.1 


2.5E-15 


25.29 


99.00 


4.5 


G159, L402,] 


M3323 


53.184639 


-27.880917 


3.471 spec. 


0.5-2.1 


12.7 


4.5E-16 


23.98 


23.84 


4.0 


G21,L369,XE 


M4302 


53.174389 


-27.867350 


3.610 spec. 


0.4-5.0 


30.6 


1.5E-15 


22.43 


22.44 


9.6 


MY30943,G2 


M4417 


53.097237 


-27.865793 


3.470 spec 


0.5-3.2 


3.1 


2.2E-17 


23.33 


23.15 


2.5 


MY31136 


M4835 


53.078468 


-27.859852 


3.660 spec. 


0.8-6.3 


12.4 


1.2E-16 


25.08 


99.00 


7.3 


G263b,L180,: 


M5390'= 


53.124371 


-27.851633 


3.700 spec. 


0.7-6.3 


18.5 


2.3E-16 


24.62 


23.44 


6.8 


G202,L265,X 


M5522 


53.157639 


-27.849911 


3.690 (2.7-4.1) 


0.3-3.0 


3.0 


2.4E-17 


25.41 


24.31 


1.2 


MY33385 


M5842 


53.070236 


-27.845531 


3.830 (2.8-4.3) 


0.7-1.8 


3.5 


3.2E-17 


25.20 


24.68 


3.5 


X235, MY34 


M7305 


53.072567 


-27.827642 


3.440 (3.1-3.7) 
3.190 (2.5-3.9)^* 


0.9-1.8 


3.5 


3.5E-17 


24.37 


24.45 


1.7 




M8039 


53.034122 


-27.817532 


0.5-5.0 


3.1 


1.7E-17 


26.95 


-26.24 


0.1 


MY38702 


M8273 


53.165268 


-27.814058 


3.064 spec. 


0.9-6.3 


31.5 


8.0E-16 


24.51 


22.40 


16.7 


MY39298,G2 


M8728 


53.069000 


-27.807203 


4.820 (3.9-5.5) 


0.8-2.2 


3.1 


5.4E-17 


26.58 


-26.24 


0.7 




M9350 


53.167900 


-27.797960 


3.360 (3.2-3.5) 
3.610 (3.0-4.1)'' 


0.4-3.6 


3.1 


1.7E-17 


24.76 


24.74 


0.8 


MY41761 


M9363 


53.051006 


-27.797874 


1.3-4.0 


3.2 


8.4E-17 


26.82 


-26.24 


0.1 




M10390 


53.059620 


-27.784552 


3.210 (2.4-3.7) 


0.3-3.6 


2.8 


2.0E-17 


26.32 


24.90 


1.0 




M10429 


53.178471 


-27.784035 


3.193 spec. 


0.4-4.5 


22.7 


7.6E-16 


25.18 


24.08 


2.0 


MY44068,L3 


M10515 


53.143574 


-27.783009 


5.260 (1.8-6.6) 


1.0-6.0 


3.3 


7.3E-17 


25.40 


23.42 


21.1 




M10548 


53.021168 


-27.782366 


4.823 spec. 


0.4-1.9 


3.2 


4.4E-17 


23.71 


99.00 


3.4 


MY44369 


M13473 


53.171535 


-27.743347 


3.210 (2.6-3.8)'' 


0.6-1.8 


3.1 


2.6E-17 


24.80 


24.23 


0.7 


MY50924 


M13549 


53.020576 


-27.742151 


3.462 spec. 


1.9-6.3 


3.7 


1.2E-16 


23.69 


23.33 


3.5 


MY51260,E1 


M14078 


53.117867 


-27.734310 


3.560 (2.8-3.8) 


0.3-4.0 


5.2 


4.7E-17 


23.97 


22.69 


19.9 


E2309,L252,: 


M70091 


53.158215 


-27.733717 


3.200 (2.0-5.4) 


0.9-2.0 


4.2 


4.2E-17 


27.65 


24.55 


5.3 


E2551,G557,: 


M70099 


53.119797 


-27.743099 


6.220 (> 2.9) 


0.5-1.9 


3.7 


3.7E-17 


27.10 


24.91 


6.2 


X392,E1516 


M70168 


53.111530 


-27.767839 


5.010 (2.0-6.5) 


0.8-4.7 


11.1 


1.3E-16 


25.97 


-25.03 


2.3 


G518,L245,X 


M70340 


53.040970 


-27.837717 


3.900 (> 3.9) 


0.4-4.2 


9.2 


l.lE-16 


-27.56 


-25.46 


2.2 


G537,L103'', 


M70385 


53.107418 


-27.855694 


3.980 (1.1-4.7) 


0.6-2.0 


7.3 


1.2E-16 


26.49 


-25.26 


4.1 


G589,L235-'',: 


M70390 


53.085411 


-27.858023 


4.180 (2.2-4.8) 


0.5-1.7 


3.8 


3.8E-17 


26.20 


23.98 


5.7 


L191,X285 


M70407 


53.064629 


-27.862434 


5.900 (>1.8) 


0.9-1.9 


2.8 


2.8E-17 


27.02 


24.89 


1.5 




M70429 


53.137917 


-27.868238 


3.630 (1.2-4.8) 


0.4-3.4 


9.7 


1.5E-16 


-27.90 


25.41 


2.3 


G217a,L290® 


M70435 


53.215118 


-27.870247 


3.230 (>2.0) 


0.7-6.0 


9.8 


1.2E-16 


26.88 


24.06 


4.9 


G508,L404,X 


M < U4o / 


0o.14d4d1 


-Z/.o / U940 


A A on ^ ^ \ 
4.42U I ) 




1 .1 


Q OTP 1 7 


/< .DZ 


-zo.oO 


0.0 


T QAa'^ W AQ^ 


M7n4fi7 




-27 882565 


4 250 ?>2 


0.9-3.8 


3.0 


1.8E-17 


26.85 


-25.50 


3.9 




M70481 


53.189320 


-27.888466 


4.080 (3.3-4.5) 


0.3-6.7 


2.9 


l.lE-16 


26.75 


24.02 


8.0 




M70508 


53.210999 


-27.907394 


3.180 (2.0-4.0) 


0.3-5.7 


4.6 


9.7E-17 


26.24 


99.00 


7.5 




M70525 


53.136761 


-27.917326 


3.280 (0.8-3.8) 


0.9-6.7 


3.1 


1.2E-16 


-27.87 


24.26 


3.7 




M70531 


53.151512 


-27.926792 


3.260 (>1.9) 


1.7-4.5 


3.0 


l.lE-16 


26.07 


23.27 


19.7 




M3607'' 


53.215858 


-27.876823 


3.468 spec 


1.0-3.8 


2.9 


8.5E-17 


24.93 


99.00 


0.9 




M70107'' 


53.079334 


-27.741474 


3.356 spec 


0.8-4.7 


3.4 


2.1E-17 


26.85 


23.77 


12.8 


E1611 


216554'''' 


53.068958 


-27.684222 


>6.0 


0.4-0.8 


2.4 


1.4E-5' 


25.8 









''L=Luo08, M=GOODS-MUSIC, MY=GOODS-MUSYC Cardamone et al. 2010;X=Xuell,G=Giacconi et al. 2002; Gilh et al. 
2011; " Norman et al. 2002; ''secondary solution at z<0.5; " LuoOS z=2.45; Luo08^ 2.71; ^ Luo08 z=2.36; Luo08 z> 7.4; just 
below detection threshold; ' Bouwens et al. 2006; ' count rate. 



A(0.16/im) lower limit obtained using the la lower limit to 
the neutral gas absorption column density and assuming a 
Galactic dust-to-gas ratio. We also plot for these sources 
the best fit N/j and A(0.16/zm) upper Hmits in Fig. [3] We 
see that measurements and the A(0.16/im) limits of all 
ten highly obscured or Compton thick AGN are completely 
inconsistent with a Galactic dust-to-gas ratio. This is true 
also for the four objects with point like morphology. The 
same conclusion is reached by estimating dust extinction 
from the galaxy stellar masses. It is difficult to estimate 
accurate stellar masses for our high-z AGN, but they likely 
exceed several 10^° Mq. The extinction expect ed based on 
the correlations found by IPannella et al. 20091 in a sample 
of z~ 2 BzK galaxies is of 3-6 magnitudes at 1500 A . We 



conclude that, at least for the four point-like objects, the 
dust-to-gas ratio of the absorbing matter is likely 1/100- 
1/1000 that of the Galaxy. This is also lower than in the 
(mostly low-z) AGN studied by IMaiolino et al. 20011 and 
IShi et al. 20061 The low dust-to-gas ratio may be explained 
if the absorbing matter is within of close to the dust subli- 
mation radius. The other six highly obscured or Compton 
thick AGN have extended morphologies in the z and H im- 
ages, suggesting that young stellar populations significantly 
contribute to the UV and optical rest frame emission. In 
these cases little can be said about the nuclear extinction 
and the dust-to-gas ratio. 
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Table 3. Multiwavelength properties of z > 3 galsLxies 



ID 


r 


logL(2-10keV) 


X-ray 


F(1.4GHz) 


SFR 


Morph.'' 


A(0.16) 


A(0.16) 


Optical 




X-ray 


ergs/s 


obsc." 


^Jy 


Mg/yr 


H band 


SED 


X-ray 


spec. 


E537 


-0.15±0.2 


43.8 


CT 


19.6±6.2 


730±230 


E 


< 2.6 


> 1000 




E737 


1.6±1.0 


43.1 




15.0±6.2 


360±150 


E 








E1312 


1.2±0.7 


43.7 








P 








E1516 




42.8 




17.8±6.4 


660±240 


E 








E1577 


-o.9tii 


43.6 


CT 






P 


< 1.1 


> 930 


Em.: Lya narrow, CIV 
Abs.: CIV broad 


E1617 


0.9±0.3 


44.4 


HO 






E 


< 2.6 


> 190 




E2199 


1.6±1.3 


42.9 








P 








E2309 


1.7±0.7 


42.9 








E 








E2498 


0.9±0.5 


43.9 




16.3±6.4 


1500±590 


E 








E2551 


1.8±0.7 


43.2 




33.5±6.49 


1870±360 


E 








E2658 


0.3±1.4 


42.8 








P 








E4956 


2.2±0.9 


43.9 








E 








E5165 


1.8±0.8 


43.2 








E 








E6257 


1.5±1.5 


42.8 








E 








E7911 


1.5±1.0 


43.6 








P 








E8479 


-0.1±0.9 


43.5 


CT 






E 


< 1.1 


> 470 




E10247 




43.1 








E 








ID 


r 


logL(2-10keV) 


X-ray 


F(1.4GHz) 


SFR 


Morph." 


A(0.16) 


A(0.16) 


Optical 




X-ray 


ergs/s 


obsc." 


/ijy 


Mo/yr 


z band 


SED 


X-ray 


spec. 


M208'^ 


0.1±0.6 


44.0 


CT 


19.4±6.5 


900±300 


P 


< 2.0 


> 780 


Em: Lya narrow, NV bro 


M3320 


1.12±0.06 


44.8 


HO 






E 


< 3.5 


> 120 


Em.: Lya narrow 


M3323 


1.6±0.21 


43.8 








P 






Em.: Lya broad 


M4302 


1.33±0.08 


44.4 


HO 


17.6±6.1 


450±160 


P 


< 1.8 


> 110 


Em.: Lya, Ly/3, Ly7, 
Abs: NV, CII, OI broad 


M4417'' 


1.8±1.0 


42.5 




15.3±6.3 


360±150 


E 






Abs: CIV, Sill, SilV, 
CII, 01, Sill 


M4835 


0.08±0.24 


44.0 


CT 


36.0 ±6.2" 


960±160 


E 






Em.: Lya, CIV narrow 


M5390'= 


0.18±0.16 


44.3 


CT 


19.4 ±6.1 


530±170 


P 


< 1.1 


> 820 


Em.: Lya, CIV, OVI, 
NV, CIV, Hell narrow 


M8273 


0.29±0.10 


44.8 


HO 


49.6 ±6.2" 


900±110 


E 


< 2.0 


> 380 


Em.: Lya, CIV narrow 


M10429 


1.40±0.10 


44.0 








P 






Em.: CIV broad 


M10548 


2.75±1.30 


43.0 








E 






Em.: Lya narrow, noisy 


M13549^ 


-0.9tr4 


43.6 


CT 






P 






Em.: Lya narrow, CIV, 
Abs: CIV broad 



" CT=Compton thick (Njj>10^'* cm"^); HO=HighIy obscured, Compton thin,(10^^ S Nh < 10^'' cm"^); " Morphology, P=point 
hke, E=Extended; " Gilh et aL (2011); Maiohno et al. 2008, [OH], [OIII] H/3 " Norman et al. (2002); ^=E1577; "VLA-CDFS 
DR2 flux 



3.2. Radio counterparts and star-formation rates 

Two Chandra-GOODS-ERS sources (E2551 and 
E1611/M70107) and other five Chandra-GOODS-MUSIC 
sources (M2690, M4835, M8273, M70091 and M70340) 
have a detection at 1.4 GHz in the DR2 catalog of the 
VLA-CDFS survey Mil ler et al.^0081 Radio fluxes for 
the Chandra-GOODS-ERS sources and the Chandra- 
GOODS-MUSIC sources with spectroscopic redshift are 
given in Table 3. Other 8 sources in Table 3 have a faint 
radio signal at the position of the galaxy (signal to noise 
between 2.4 and 3.2). The probability that this signal is 
a background fluctuations is < 2%, corresponding to < 1 
spurious radio detection in the full sample of sources in 
Table 3. We report in Table 3 radio fluxes also for the 
additional 8 faint detections. We stress that the latter 
faint fluxes may over-estimate the real radio flux because 
of the Eddington bias. A more robust average flux can 
be obtained by stacking together the radio images at 
the position of the X-ray sources. The average flux for 
the Chandra-GOODS-ERS and Chandra-GOODS-MUSIC 
sources without a detection in the DR2 VLA-CDFS catalog 
is 6.5±1.5/xJy and 10.1±l.l/iJy respectively. 



The observed radio fluxes can be due to both nuclear 
and stellar processes. Using the lLa Franca et al. 201l |prob- 
ability distributions for the X-ray to radio nuclear lumi- 
nosity we would expect a L(2-10keV)/L(1.4 GHz) ratio 
< 3 X 10"* for ~ 25% of the sources in Table 3, and L(2- 
10keV)/L(1.4 GHz)< 10^ for 5% of the sources. We con- 
verted the radio fluxes in Table 3 into luminosities by using 
the following expression: 

L(1.4GHz) ^ 1.19 X 102"i:)L2/(1.4GHz)(l ± z)("~i) (1) 

where DL is the luminosity distance in Mpc and 
/(1.4GHz) is in /iJy and we used a — 0.7. We flnd in all 
cases L(2-10keV)/L(1.4 GHz)< 3 x IQ-*, which is 41% of 
the fufl sample, and L(2-10keV)/L(1.4 GHz)< 10^ for 5 
sources, i.e. 18% of the full sample. This difference between 
the expected and observed fraction of radio bright sources 
is not conclusive, because our radio fluxes may be over- 
estimated as discussed above. However, if we associate to 
the sources without a detection in the DR2 VLA-CDFS 
catalog the average flux measured above, the fractions of 
sources with L(2-10keV)/L(1.4 GHz)< 3 x lO'' and L(2- 
10keV)/L(1.4 GHz)< 10^ are 92% and 19% respectively. 
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1 2 

A(0.16/xm) 

Fig. 3. The best fit N/f as a function of the 0.16/ini ex- 
tinction for the 10 z> 3 highly obscured or Compton thick 
AGN. Filled points are sources with point like optical or 
near infrared morphology. The two solid lines are the ex- 
pectation for a dust-to-gas ratio 100 and 1000 times lower 
than the Galactic one. 

We can therefore conclude that, at least on average, radio 
fluxes are higher than expected based on the X-ray fluxes 
and assuming common nuclear origin. If radio emission is 
dominated by stellar processes we can convert radio fluxes 
in star- formation rates (SFRs, see e.g. lYun et al. 200T|) . We 
used the following calibration which assumes a Chabrier 
IMF0: 



SFR = 3.4 X IQ-^^LilAGUz) 



(2) 



In two cases we have independent estimates of SFR, 
which agrees reasonably well with the values in Table 3. 
M208, has a SFR~ l OOOMg/yr estima ted from a Laboca 
870^m detection (see IGilh et al. 201 1|) . while M4417 has 
SFR- 44OM0/y r, estimated fro m UV SED fitting and oxy- 
gen lines (|Maiolino et al. 2008p . 

3.2.1. A high-z star-forming galaxy? 

Rather peculiar is the case of M4417, a galaxy studied 
by IMaiolino et al. 20081 in the framework of the AMAZE 
program. The X-ray luminosity corresponding to the 
observe d SFR is logL(2-10keV)- 42.35 (by using the 
[Ranalli et al. 20031 conversion factor), just a factor 40% 
lower than the observed X-ray luminosity. It should be 
noted that although the Chandra X-ray contours seems 
centered on M4417 (see Fig. 2]), there could be some contri- 
bution to the observed X-ray flux from the nearby z=3.471 
galaxy M4414, which has a SFR - 1/4 that of M4417. 
Intriguingly, M4417 is the Chandra-GOODS-MUSIC and 
Chandra-GOODS-ERS object with the lowest X-ray to z 
band and H band flux ratios (0.004 and 0.005 respectively), 
lower than typical AGN (see next section) , further suggest- 
ing that the X-ray emission from this object is not dom- 
inated by a nuclear source. We then conclude that the 



observed X-ray luminosity is likely dominated by stellar 
sources, making M4417 one of the farthest objects in which 
X-rays are probing stellar processes. Of course other solu- 
tions are possible, for example M4417 may be a reflection 
dominated, heavily Compton thick AGN, where the nucleus 
is completely hidden and the observed X-ray luminosity is 
just a fraction of the real one. We do not have however ev- 
idences in this direction from either X-ray colors, optical 
spectroscopy, infrared spectroscopy and broad band UV to 
infrared SED (the source is not detected at wavelengths 
longer than 8/xm, see Fig. S]). 

There are other two sources with relatively low X-ray 
luminosity and high radio flux: E2551 and E1516. In these 
cases the X-ray luminosity corresponding to the SFR rate 
estimated from the radio flux is a factor of 70% and 90% 
lower than the observed X-ray luminosity. Their X-ray to 
H band ratio are however Z 10 times higher than that of 
M4417 (0.06 for E2551 and 0.07 for E1516). The radio flux 
and luminosity of E2551 are the highest in the sample, sug- 
gesting that this may well be a radio loud AGN. 



3.3. 



X-ray to optical/NIR flux ratios and selection effects 



^ assuming a Salpeter IMF would result in SFR 
higher 



1.7 times 



The standard procedure to assemble X-ray high-z AGN 
samples is through the optical identification of X-ray source 
catalogs. We used in the previous section a complementary 
approach, that is studying the X-ray emission of optically or 
NIR selected high-z galaxies. In the first case the problem is 
to be sure that the optical/NIR identification of the X-ray 
source is correct, and that sample identifications are rea- 
sonably complete. In the latter case the problem is the reli- 
ability and completeness of the optical/NIR samples. They 
should include most of the galaxies (and therefore most 
AGN) down to a reasonably low fiux limit. The GOODS- 
MUSIC catalog is fairly complete down to z '~ 26 and 
F(4.5/im)— 1.5/iJy and includes sources down to z ~ 27 
and F{4.5^iTn)^ 0.5^i3y. The GOODS-ERS catalog is fairly 
complete down to H— 26.5 and includes sources down to 
H'^ 27.5. To understand how these optical and NIR flux 
limits translate into X-ray flux limits we compute the X-ray 
(0.5-2keV) to NIR (H band) fiux ratio F(0.5-2keV)/F(H) of 
much brighter, and therefore likely complete, AGN samples. 
Fig. [5] compares F(0.5-2keV)/F(H) of the z> 3 Chandra- 
GOODS-ERS and Chandra-GO ODS-MUSIC AGN sam- 
ples to the Cha ndra-COSMOS (fCivano et al. 201ip and 
XMM-COSMOS (IBrusa et al. 2009allB?^ et al. 2010|1 z> 
3 AGN samples. The COSMOS samples are X-ray se- 
lected samples with nearly complete optical identification, 
thanks to the fact that they cover X-ray and optical/NIR 
fiux ranges 10-100 times brighter than the Chandra- 
GOODS-ERS and Chandra-GOODS-MUSIC AGN sam- 
ples, and the massive photometric campaigns performed 
in this field with HST/ACS, Subaru and CFHT (Capak 
et al. 2008). As an example, there are only two sources 
in the Chandra-COSMOS catalog without an optical 
or infrared counterpart, 19 without an optical counter- 
part, 0.1% and 1% of aU Ch andra-COSMOS sources re- 
spectively, ICivano et al. 20111 Most of these sources are 
likely high luminosity highly obscured type 2 QSOs (e.g. 
IFiore et al. 20 03). But even assuming that all these 19 
sources are at high-z, they would be 25% of the z> 3 sources 
in Chandra-COSMOS. So this is a very conservative up- 
per limit to the Chandra-COSMOS completeness at high-z. 
Similar numbers apply to the XMM-COSMOS survey. 
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Fig. 4. K band, MIPS 24/iin and LABOCA 870/zm images around the M4417 galaxy with overlaid X-ray 0.5-3.2 keV 
contours. 



The Chandra-GOODS-ERS detections cover a range of 
F(0.5-2keV)/F(H) similar to the COSMOS samples and 
indeed the Chandra-GOODS-ERS and COSMOS distri- 
butions of F(0.5-2keV)/F(H) are perfectly consistent (the 
probability that the Chandra-GOODS-ERS and Chandra- 
COSMOS distributions are drawn from the same parent 
population is 50%, using the Kolmogorov-Smirnov test). 
This suggest that the HST/WFC3 ERS H band images are 
deep enough to trace high-z AGN populations at the ex- 
tremely faint X-ray flux limits reached by the Chandra 4 
Msec exposure, avoiding significant incompleteness, or at 
least with an incompleteness comparable to that reached 
at much brighter fluxes by Chandra-COSMOS and XMM- 
COSMOS. The same exercise performed using the GOODS- 
MUSIC z band and IRAC A.5^m selected catalog is not 
satisfactorily (we used to this purpose a subsample of the 
full Chandra-GOODS-MUSIC sample with deep H band 
VLT/ISAAC coverage, excluding 13 sources with shal- 
lower H band coverage). The Chandra-GOODS-MUSIC 
sample clearly misses sources with high F(0.5-2keV)/F(II) 
flux ratio in comparison to the Chandra-GOODS-ERS 
and COSMOS samples, and indeed the Chandra-GOODS- 
MUSIC F(0.5-2keV)/F(H) distribution is inconsistent with 
the Chandra-COSMOS distribution at the 99.925 % con- 
fidence level (using the Kolmogorov Smirnov test). This 
means that the GOODS ACS and IRAC images are not 
deep enough to fully probe the AGN population at the 
fiux limits of the CDFS 4 Msec observation. The F(0.5- 
2keV)/F(H) distribution of the Chandra-GOODS-MUSIC 
sources with F(0.5-2keV) 5 x 10"^^ erg/cm^/s is consis- 
tent with the GOODS-ERS and COSMOS ones (proba- 
bility of ^ 30% that they can be drawn from the same 
parent population) and therefore in the following we will 
consider only the 23 Chandra-GOODS-MUSIC sources 
brighter than this flux limit (and outside the ERS area). 
These 23 sources are added to the 17 sources of Chandra- 
GOODS-ERS sample to form a total sample of 40 high-z 
AGN. We use this sample to constrain the faint end of the 
AGN luminosity function at high redshift. 

4. High-z AGN luminosity functions 

Our z> 3 AGN candidate sample includes 17 Chandra- 
GOODS-ERS sources and 23 Chandra-GOODS-MUSIC 
sources spanning a range in luminosity 42.5<logL(2- 



10keV)/ergs/s<44.8. They can be used to probe the faint 
end of the high-z AGN luminosity functions. We calculated 
the comoving space densities of our high-z sample using 
the 1/Vmax method (jSchmidt 1968) . We choose L(2-10keV) 
and z ranges to ensure completeness at the X-ray flux limit 
F(0.5-2keV)~ 2 x 10~^^ erg/cm^/s reached by our survey. 
The luminosity limit at z=4, 5 and 7.5 is logL(2-10keV)'-~^ 
42.7, 42.8, 43.3 respectively. Accordingly we compute co- 
moving space density in the redshift bins 3-4, 4-5 and 5.8- 
7.5 with luminosity ranges 42.75-44.5, 43-44 and 43.5-44.5 
respectively. There are a total of 30 Chandra-GOODS- 
ERS and Chandra-GOODS-MUSIC AGN in these redshift 
and luminosity bins. Comoving space densities are given in 
Table 4. Fig. [7] presents the AGN luminosity functions in 
the same redshift bins. Errors are computed using Poisson 
statistics on the number of AGN in each redshift-luminosity 
bin. To obtain information on the slope of the luminosity 
functions we joined the above samples with the Chandra- 
COSMOS high-z sa mple of[C.i vano et al. 2"0Tn the XMM- 
COSMO S sample of [Bmii| iral. 2009a, the GOODS sam- 
ple of Fonta not et al. 20071 the faint optical AGN sample 
of Glikma n et al. 20111 the luminous optical AGN sam- 
ples of lRichards et al. 20061 and [Jiang et al. 2009] We con- 
verted the rest frame 1450 A luminosities of the opti- 
cally selected samp les to the 2-10 keV band using the 
M arconi et al. 2004l and Sirigu et al. 2011 luminosity depen- 
dent conversion factors. We assume no intrinsic reddening 
in the optically selected, high-z AGN. This might underesti- 
mate the real AGN luminosity, if signiflcant dust is present 
along the high-z AGN lines of sight (M aiolino et al. 20041 
Jiang et al. 2006[ IGallerani et al. 2010|l . We checked our 
conversion factors against real data, using the Chandra 
and XMM detecti o ns of the [Jiang et al. 2009| QSOs 
([Mathur et al. 20"02l IShemmer et al. 20061 and references 
therein). Fig. [HI shows the ratio between the UV 1450 A 
luminosity and the 2-10 keV luminosity against 1450A lu- 
minosity for 12 QSOs at z> 5.7 for which we collected X-ray 
data from the literature. The flgure also shows the UV to 
X-ray conversion adopted in this work, which agrees quite 
well with the data of this high redshift QSO sample. 
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Fig. 5. [Left panel]: The X-ray to NIR flux ratio F(0.5-2keV)/F(H) as a function of the X-ray flux for the z> 3 Chandra- 
GOODS-ERS AGN sample (red-filled circles), Chandra-GOODS-MUSIC sample (black-open circles), Chandra-COSMOS 
sample (green-filled squares) and XMM-COSMOS sample (blue-filled triangles. [Right panel:] The distributions of F(0.5- 
2keV)/F(H) for the Chandra-COSMOS z> 3 sample (black histogram), Chandra-ERSVO sample (red histogram), and 
Chandra-GOODS-MUSIC sample (blue, dashed histogram). 



Table 4. Chandra-GOODS-ERS z>3 AGN comoving space densities 
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Fig. 6. The UV 1450A to X-ray 2-10 keV luminosity ratio 
for 12 QSOs at z> 5.7 

We modeled the wide luminosity range high-redshift lu- 
minosity functions using the standard double power law 
shape: 

d3'(ix) 



where A is the normalization factor for the AGN den- 
sity, 7I and 72 are the faint-end and bright-end slopes and 
L* is the characteristic break luminosity. Optical selection 
may miss highly obscured AGN, which can represent a 
large fraction of the total in particular at low luminosity 
(|La Franca et al. 2005p . To avoid possible incompleteness 
in optical surveys we excluded from the fit the optically 
selected AGN density determinations with M1450 > —26.5 
(or correspondingly L(2-10keV)<45.1). We first fitted this 
simple model to the luminosity functions in the three red- 
shift bins. The number of points in the z=3-4 redshift bin 
allows us to constrain all four parameters. This is not pos- 
sible in the redshift bins 4-5 and > 5.8, where we fitted the 
data fixing the slopes 7I and 72 to the best fit values found 
in the z—S-A redshift bin. The results of these fits are re- 
ported in Table 5. As a next step, we fitted simultaneously 
the data in the three redshift bins with an evolutionary 
model. We chose to describe the evolution of the high-z 
AGN luminosity function with the LADE (luminosity and 
density evolution) model, introduced by lAird et al. 20101 
Thus the evolution of L* is given by: 



LogL,(z) = LogLo +pLog(l + z) 

and the evolution of A is given by: 
LogA(z) = Log^o -I- c?Log(l + z) 



(4) 



(5) 



dLogL> 



A[{L^/L,y^ + (Lx/L,)^'] 



(3) 
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Fig. 7. The AGN 2-10 keV luminosity functions at 
z=3-4 (top panel), z=4-5 (central panel) and z>5.8 
(bottom panel). Large red-filled circles are from the 
CDFS Chandra-GOODS-MUSIC and Chandra-GOODS- 
ERS samples. Red-filled triangles are from Chandra- 
COSMOS (Civano et al. 2011), red-filled squares from 
XMM-COSMOS (Brusa et al. 2010). The blue thick curve 
in the top panel is the LADE model of Aird et al. 2010 in the 
2.5-3.5 redshift bin. Large black-open triangles are SDSS 
data from Richards et al. 2006, stars are SDSS data from 
Jiang et al. 2009, open circles are from Glikman et al. 2011 
and small blue-open triangles are from Fontanot et al. 2007. 
Black dashed curves are best fits to the X-ray plus optical, 
wide luminosity range luminosity functions (see text for de- 
tails). Magenta solid curves encompass predictions on the 
Menci et al. (2006, 2008) semi-analytic model. Blue curves 
are models from Shankar et al. (2011), see the Discussion 
for details. 



We first fitted the data with the full six parameter 
model. The 7I and 72 slopes turned out similar to those 
obtained fitting the z=3-4 data only (the z=3-4 data pro- 
vides the strongest constrain to the shape of the luminosity 
function). We report in Table 5 the best fit parameters ob- 
tained by fixing 72 to 3.4 and to 3.0 (keeping the faint-end 
slope 7I fixed at 0.8 in both cases). The fit with 72 = 3 
produces a significantly worse than for the best fit case, 
but still acceptable considering the uncertainties on the cor- 
rection between the optical and X-ray luminosities for the 
optically selected AGN density determinations. We found 
that the Aq and d parameter are completely degenerate. 
We could obtain equally good fits with combinations of Aq 
and d inversely correlated. We report in Table 5 the best 
fits obtained with d = (no density evolution) and d = —1. 
On the other hand, the data are good enough to constrain 
relatively well Lq and p. Our data are thus consistent with 
a pure luminosity evolution, with rather quickly reduc- 
ing with the redshift. The dashed lines on Fig. [7] are the 
best fit pure luminosity evolution model (4) in Table 5. 

It is not straightforward to compare our high-z, wide lu- 
minosity band luminosity functions with previous determi- 
nations. The best fit model of lShankar et al. 2009al agrees 
quite well with our determinations at logL(2-10keV)< 44 
but overestimate the density of higher luminosity AGN, in 
particular of optically selected AGN. It should be noted 
that that IShankar et al. 2 009a adopt UV to X-ray lumi- 
nosity conversion factor fixed at 10.4, while our conver- 
sion factor varies with the luminosity, from 4.4 at logL(2- 
10keV)=43 to 27 at logL(2-10keV)=45. Furthermore, 
IShankar et al. 2009al correct for extinction and the frac- 
tion of Compton thick AGN. As discussed above our data 
are not corrected for the fraction of Compton thick ob- 
ject, because several have actually been found in our sam- 
ples, and we use optically selected AGN of high luminos- 
ity (M1450 > -26.5, L(2-10keV)<45.1) where the fraction 
of obscured AGN is likely to be small. [Hopkins et al. 2007| 
use a luminosity dependent bolometric correction, but with 
a different calibration with respect to ours. They correct 
their data for extinction and the fraction of Compton 
thick AGN missed in optical. X-ray and infrared surveys. 
They convert Njj distributions evaluated from X-ray data 
into optical-UV reddening using a canonical gas-to-dust 
ratio, which, as discussed in the previous section (and 
by e.g. IMaiolino et al. 20011 IShi et al. 20"06|l . can overes- 
timate the real extinction, in particular when the X-ray 
absorber is compact, smaller than the dust sublimation ra- 
dius. They find z> 3 luminosity functions with much flatter 
slopes with respect to ours (or to the [Shankar et al. 2009al 
ones). As an example, the [Hopkins et al. 2007| best fit 
model would predict ~ 2 — 4 times more z=3-4 and z=4-5 
AGN with logL(2-10keV)> 44.5 than a ctually found in the 
XMM and Chandr a COSMOS surveys (jBrusa et al. 2009al 
ICivano et al. 201 Ij) . 

4.1. AGN luminosity function and duty cycle evolution 

It is interesting to put the findings of the previous section 
in a context, to study the evolution of the AGN luminosity 
function from the local Universe to 6. We compare in 
the upper panel of Fig. [S|the best fits to the z=3.5, z=4.5 
and z=6 AGN luminosity functions (model 4 in Table 5) to 
the best fit model luminosity functions found at lower red- 
shift bv ILa Franca et al. 20051 These have been computed 
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Table 5. Modeling the high redshift AGN luminosity function 



model 


z range 


A or Ao 
10"" Mpc"^ 


I/* or I/O 
10** ergs/s 


7I 


72 


P 


d 


X" (dof) 


1 


3-4 




5.4 ±2.2 


0.8 ±0.2 








4.63 (4) 


2 


4-5 


3.2+11 


3.5 ±0.4 


0.8 FIX 


3.4 FIX 






2.52 (2) 


3 


> 5.8 






0.8 FIX 


3.4 FIX 






1.23 (1) 


4 


> 3 


q q + U.V 

'J-'J-l.O 




0.8 FIX 


3.4 FIX 


-2.1 ±0.2 


FIX 


10.65 (14) 


5 


> 3 


^■^-1.6 


1301™ 


0.8 FIX 


3.0 FIX 


-2.4 ±0.2 


FIX 


16.09 (14) 


6 


> 3 


l-"-0.5 


7011° 


0.8 FIX 


3.4 FIX 


-1.8 ±0.2 


-1 FIX 


10.82 (14) 


7 


> 3 


^•0-0.8 


751^5 


0.8 FIX 


3.0 FIX 


-2.03 ±0.15 


-1 FIX 


15.77 (14) 



by correcting for the incompleteness due to X-ray absorp- 
tion, which is important at low redshift where the photo- 
electric cut-off produced by the typical column densities 
observed in many AGN is well within the X-ray selection 
band. The iLa Franca et al. 20051 luminosity functions are 
consistent with the new determinations of lAird et al. 20101 
and lEbrero et al. 20091 based on a larger number of ob- 
jects. The upper panel of Fig. [5] summarizes the complex 
evolution of the AGN luminosity function, which rises from 
z=0 up to z=l, 2, 3 for increasing luminosities and than 
decreases at higher redshifts. 

We converted these luminosity functions into "active" 
SMBH mass functions by using MonteCarlo realizations. 
We simulated 10^ AGN luminosities and SMBH masses in 
each redshift bin according to the following procedure. 1) 
We first randomly choose an X-ray luminosity following 
the luminosity function distribution in each given redshift 
bin; 2) we than convert it into a bolometric luminosity us- 
ing the IMarconi et al. 20041 and Sirigu et al. 2011 lumi- 
nosity dependent bolometric correction; 3) we randomly 
choose an Eddington ratio from given distributions. We 
used log-normal distributions with parameters in Table 
6. At z< 0.3 we used the distribution of INetzer 2009b[ 
which are shifted toward higher Eddington ratios with 
respect to the IKauffmann fc Heckman 2009' distribu- 
tions (see the discussion in jNctzcr 2009b, Nctzcr 2009a|. 
At medi um to high redshift w e used the distribu- 
tions of jTr akhten brot et al. 20111 IShemmer et al. 20041 
INetzer fc Tfa khtcnbi^ oFSOOTl IWillott et al. 2010bl This 
provides a SMBH mass for each chosen X-ray luminos- 
ity and Eddington ratio. 4) We finally binned the re- 
sulting SMBH distributions in each redshift bin to build 
SMBH mass functions. The resulting "active" SMBH mass 
functions in six redshift bins are plotted in the central 
panel of Fig |S| It must be stressed that this is an empir- 
ical calculation, performed independently in each redshift 
bins, using only observed quantities. It does not pretend 
to model the history of SMBH growth through the cos- 
mic times, which can be obtained by using a continuity 
equation, conserving the number, as in lMarconi et al. 20041 
iMcrloni fc Heinz 2 008. Shankar et al. 2009al 

Furthermore, the adopted Eddington ratio distributions 
all refer to relatively luminous AGN. Less active SMBH, 
giving rise to low luminosity AGN and LINERs are not rep- 
resented in these distributions. For this reason we adopted 
a luminosity limit for the computed SMBH mass func- 
tions: these functions represent "active" SMBH produc- 
ing an X-ray luminosity > 10'*'^ ergs/s. As a consequence, 
comparison with previous calculation is not straightfor- 
ward. For example, while our SMBH mass functions are 
luminosity limited, the "active" SMBH mass functions of 



'Me rloni fc Heinz 20081 are X-ray flux limited mass func- 
tions. The 2-10 luminosity corresponding to the faintest flux 
limit of IMerloni fc Heinz 20081 is ~ 10''° ergs/s at z=0.1, 
- 10^^-^ and ~ lO'^^-S at z=3, while we plot the SMBH mass 
function of AGN more luminous than ^ 10*'^ in all redshift 
bins. SMBH mass functions of broad line AGN have been 
computed by Kelly et al. 2010, by using about 10.000 SDSS 
QSOs in the redshift range 1-4.5. Our SMBH mass func- 
tions are consistent, or slightly higher than the Kelly et 
al. determinations at most redshifts and black hole masses 
considered, as expected since broad lines AGN are a frac- 
tion of the full active SMBH population, even at the highest 
SMBH masses and/or AGN luminosities. The only regime 
where our estimates are signiflcantly higher than the Kelly 
et al. ones are the highest masses (logM> 9.3M0) at z~ 1, 
where the Kelly et al. 2010 function drops steeply, while 
our function decreases more smoothly. 

The SMBH mass functions can be transformed into stel- 
lar mass function of "active" galaxies by assuming a con- 
version factor between SMBH mass and host galaxy stellar 
mass. We assu med a mean value F q = log{MBH/M^) = 
-2.8 at z ^ (jHaring fc Rix 2004| but also see the dis- 
cussion in ILamastr a et al. 2010p and a redshift evolution 
rox(l+z)°-^ (Merlon TeFal. 20041 [Hopkins et al. 2006 



IShankar 2009|) . It should be noted that in the local Universe 
the above correlation has been found and calibrated for 
the bulge component of galaxies. Whether a strict distinc- 
tion between bulge and disk still exist at high-z is a mat- 
ter of debate. Disks of z^ 2 galaxies are much more com- 
pact and much thicker that in today spirals of similar mass 
fG enzel et al. 20061 Ivan der Wei et al. 2011i . For the sake 
of simplicity we assume in the following calculation that 
the SMBH mass is proportional to half of the total stel- 
lar mass of high-z galaxies. It is instructive to compare 
the SMBH mass functions and the "active" galaxies stellar 
mass functions to the stellar mass functions of all galax- 
ies. The lower panel of Fig. [S] plots a collections of galaxy 
stellar mass functions: we used the stellar mass functions 
oflFontana et al. 2006l at z < 2, ISantini et al. 2011la t z=3-4, 
jCaputi et al. 2011| at z=4-5 and 'S tark et al. 2009 ' at z=6. 

The AGN fraction, or AGN duty cycle, can finally be 
obtained by dividing the "active" galaxy stellar mass func- 
tions by the galaxy stellar mass functions. Following the 
adopted luminosity limit used to compute SMBH mass 
functions, we define AGN duty cycle the fraction of AGN 
with 2-10 keV luminosity greater than 10^'^ ergs/s to the 
total number of galaxies with a given stellar mass. Lower 
luminosity AGN do exist and are probed by the luminos- 
ity functions in Fig. |8| up to z~ 2 — 3. However, they are 
below the fiux limit of current surveys at higher redshift. 
Therefore, our luminosity threshold also avoids problems of 
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Table 6. Parameters of log-normal Eddington ratio 
distributions 



0.001 



z 


peak 


width 






dex 


0-0.3 


0.03 


0.5 


0.3-1 


0.1 


0.33 


1.5-2.5 


0.18 


0.33 


3-4 


0.22 


0.33 


4-5 


0.5 


0.33 


> 5.8 


1.0 


0.25 



incompleteness of the samples at high-z. According to this 
definition, the AGN duty cycle can be different from the 
AGN timescale, since the latter is correlated with the total 
inntrinsic lifetime of the AGN, which can include phases 
with luminosity < 10*'^ ergs/s. The AGN duty cycle is 
plotted in the upper panel of Fig. |9l It must be stressed 
again that this is an empirical calculation, performed in- 
dependently in each redshift bin, using only the observed 
AGN and galaxy luminosity and mass functions and ob- 
served Eddington ratio distributions. We find that the AGN 
duty cycle increases at all redshift with the stellar mass. 
Similar trends have been found bv IKauffmann et al. 20031 
[Best ct al. 2005', |Bundy et al. 20081 lYamada et al. 20091 
[Brusa et al. 2009b|for optically selected, radio selected and 
X-ray selected AGN. At each given stellar mass the duty 
cycle increases with the redshift up to z=4-5, consistent 
with the result of iMarconi et al. 20041 IBrusa et al. 2009b[ 
IShanka r et al. 2009al There are 22 galaxies in the GOODS- 
ERS catalog with z> 3 and stellar mass higher than 10^^'^^ 
Mq, 7 of which are X-ray sources in Table 1, thus con- 
firming an AGN duty cycle ^ 30% at z>3 for such massive 
galaxies. Fig. [HI lower panel, plots the evolution of the AGN 
duty cycle as a function of the redshift for two galaxy stellar 
masses: log(Ms)=11.25 and 11.75. Bands are plotted rather 
than curves to emphasize the rather large uncertainties, es- 
pecially at high-z (see next section). The same Fig. shows 
previous evaluations by IBrusa et al. 2009bl These are con- 
sistent with the present estimates within their rather large 
error bars. The expectations of the lMenci et al. 20081 semi- 
analytic model are also showed in the lower panel of Fig. 

m 



4.1.1. Error budgets 

The determination of the evolution of the AGN duty cycle is 
plagued by large uncertainties, especially at high redshifts. 
It is therefore important to study in detail the origin of 
these uncertainties and the way of reducing them. 

The uncertainty on both the AGN comoving densi- 
ties at z<2 (at least for unobscured and moderately ob- 
scured AGN) and on the stellar mass functions at z< 2 
are relatively small, of the order of 10-20% or even smaller, 
thanks to large AGN and galaxy samples used for these 
determinations and thanks to the use of several differ- 
ent surveys, which helps in reducing the systematic er- 
ror. The largest uncertainty in the AGN luminosity func- 
tion is on the fraction of Compton thick AGN. The frac- 
tion of these objects in the local Universe is high: 30- 
50% of the optically selected Seyfert 2 galaxies can be 
Com pton thick 1/4 - 1/3 of the fuU A GN popula- 
tion, IRisaliti et al. 19991 IPanessa et al. 20061 a result con- 
firmed by hard X-ray selection, see IMalizia et al. 2009|) . 
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Fig. 8. [Top panel:] A collection of AGN luminosity func- 
tions. z<2 from La Franca et al. 2005, z> 3 this work. 
[Central panel:] SMBH mass functions obtained by com- 
bining the AGN luminosity functions in the top panel with 
Eddington ratio distributions (see the text for details). 
[Bottom panel:] A collection of galaxy stellar mass func- 
tions. z< 2 from Fontana et al. 2006, z=3-4 from Santini et 
al. 2011, z=4-5 from Caputi et al. 2011, z=6 from Stark et 
al. 2009. 



At higher redshift the fraction of Compton thick AGN 
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is more uncertain but it can be as high as 1/3 of the 
full AGN population (Fiore et al. 2008', 'Fiorc ct al. 2009), 
or even higher (Daddi et al. 2007, Treister et al. 2010). A 
fraction of Conipton thick sources ~ 1/3 of the total is ac- 
tually included in the model of ILa Franca et al. 20051 and 
therefore we are confident that the error on the total AGN 
comoving space density at z< 2 due to undetected, and un- 
accounted for, Conipton thick sources is small ( < 10 — 20%) 
At z> 3 the uncertainties on both AGN luminosity 
function and galaxy mass functions are larger, in particular 
at low AGN luminosities and low galaxy masses (see Table 
3), they are of the order of 30-50% at z=3-5. At z~ 6 
the uncertainties on the faint end of the AGN luminosity 
function and on the galaxy stellar mass function are 
extremely large, a factor 100%. These uncertainties are 
dominated by statistics in the case of AGN (only few 
detections at z>4-5) and by systematics in the case of the 
galaxy mass function (at z~ 6 stellar masses are estimated 
by IStark et al. 20091 by converting the UV luminosity 
into stellar mass, which is a highly uncertain procedure). 
The uncertainty due to undetected or unaccounted for 
Compton thick AGN at z> 3 is smaller than the above 
features. In fact, we note that at z> 3 the cut-off produced 
by a column density lO^'' cm~^ is shifted below 2.5 keV, 
where the effective area of Chandra and XMM peaks. 
This, together with the extremely deep exposures of 
the CDFS, helps in detecting directly at least mildly 
Compton thick sources ( |Georgantopoulos et al. 2009^ 
Georgantopoulos et al. 20l"H Feruglio et al. 2011 



Gilh et al. 201 1[ IComastri et al. 20Tl| ; Indeed, about 
one fifth of the Chandra-GOODS-ERS sources can be 
Compton thick, as well as one fourth of the Chandra- 
GOODS-MUSIC sources with optical spectroscopy (see 
Sect. 3.2), and many are directly detected in our z>3 
search. We are therefore confident that the uncertainty 
on the total AGN comoving density linked to Compton 
thick AGN is at most 10 - 20% also at z> 3. The high 
luminosity end of our z>3 luminosity functions is probed 
by optical surveys. Obscured AGN may well be missed 
in these surveys, and therefore these determinations may 
be regarded as lower limits. We note, however, that the 
fraction of obscured AGN decreases strongly with the 
luminosity (at least at z< 2 — 3 where large area X-ray 
surveys provided samples of high luminosity QSOs), and 
therefore the error on the comoving space density of high 
luminosity AGN at z> 3 is also likely to be small (< 20%). 
We should consider the uncertainty in the conversion from 
UV to X-ray luminosities for the optically selected AGN 
luminosity functions. This uncertainty is probably of the 
order of 30% (see the discussion in Shankar et al. 2009a). 

Another relatively large source of error in the evalua- 
tion of the AGN duty cycle is represented by the uncer- 
tainty on the Eddington ratio distributions. We performed 
several tests using parameters slightly different than those 
in Table 4. For these AGN we found that the duty cycle 
(at logM,=ll) changes by a factor of 15-30% and ~ 100% 
for distribution peak and width differing by 30% from the 
assumed ones. The biggest effect is introduced by the width 
of the Eddington ratio distributions: the broader the distri- 
bution the higher is the resulting AGN duty cycle. The un- 
certainty on the Eddington ratios is particularly severe for 
low luminosity, logL(2-10keV)< 43 AGN. Since this is also, 
roughly, the luminosity limit of our z> 3 sample, we decided 
to limit the analysis to AGN with logL(2-10keV)> 43 ergs/s 



only, and assume conservatively a 50% (total) relative er- 
ror on the AGN duty cycle arising from the uncertainty on 
Eddington ratios. 

Finally, we adopted given normalization and evolu- 
tion of the SMBH to galaxy stellar mass ratio F and as- 
sumed that is is the same for each object at each given 
redshift. However, we know that, at least at low red- 
shift, the SMBH-bulge mass relation is unlikely to be uni- 
versal (jMathur et al. 20011 IMathur et al. 20111 and refer- 
ences therein). Furthermore. IBatcheldor 20101 suggest that 
selection effects are important in shaping the SMBH- 
bulge mass relation. Thus, the deviations in the SMBH 
mass-bulge mass relations resulting from galaxy mor- 
phology, orientation and selection effects can affect the 
duty cycle calculations, in particular at small masses. 
We assumed a redshift evolution (1 -I- z)^-^, consistent 
with our prese nt knowledge and with the expectation of 
several models ILamastra et al. 20101 [Hopkins et al. 2006 



IShankar et al. 2009al However, the uncertainty on this calf 
bration increases with the redshift, and can be rather large 
at z> 3 — 4. To quantify all these effects is not an easy 
task. As an example, if F is half of what assumed above, 
the duty cycle is reduced by a factor between 20% and 
100%, depending on the galaxy stellar mass and redshift. 
In summary, the bands plotted in Fig. [S] roughly account 
for the typical errors given above at each given redshifts 
and galaxy stellar masses. 



5. Discussion 

We evaluated the comoving space density of the z> 3 faint 
X-ray sources in three redshift bins: 3-4, 4-5 and > 5.8. 
The number of AGN in the three redshift bins is small, 
19, 9 and 2 respectively, and therefore we were forced to 
use relatively wide luminosity bins, to keep the statistical 
error reasonably small. In particular, the comoving space 
density at z> 5.8 and logL(2-10keV)=43.5-44.5 has been 
computed using only two sources. We stress that in one 
case a secondary solution in the photometric redshift does 
exist at Z'^ 2, so our determination is probably an upper 
limit to the true space density of low luminosity AGN at 
> 5.8. This confirms that the slope of the faint end of the 
z> 5.8 AGN luminosity function is significantly flatter than 
the bright end slope (see e.g. IShankar fc Mathur 2007^ . 
We also note that the other claimed z> 7 AGN in the 
CDFS (L306, M70437), has a 0.5-2 keV flux below our 
threshold for Chandra-GOODS-MUSIC sources, and that 
the photometric solution of GOODS-MUSIC is broader 
than lLuo et al. 20101 one, with a lower limit at z=2.7. This 
source is therefore not part of the > 5.8 sample used 
in Fig. [T] Recently ITreister et al. 20Tl1 evaluated the inte- 
grated AGN emissivity at 6 by stacking together the 
CDFS and CDFN X-ray data at the position of the z~ 6 
IBouwens et al. 20061 galaxy candidates. They obtained a 
luminosity density of 1.6 x 10"*^ ergs/s/deg^ in the 2-10 
keV band. If we integrate our 6 best fit model above 
logLx=42 we obtain a total luminosity density of 7.5 x 10'^* 
ergs/s/Mpc^ or 5.6 x lO''^ ergs /s/deg^, a value ~ 2 .8 times 
smaller than that reported by ITreister et al. 20111 It must 
be noted that our determination is based on a fit to the 
broad luminosity range luminosity function obtained by 
joining single X-ray detections in the CDFS to the UV lu- 
minosity functions, while the ITreister et al. 20TT] is based 
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Fig. 9. [Top panel:] the AGN duty cycle as a function of the 
galaxy stellar mass in six redshift bins. AGN with logL(2- 
10keV)> 43 are considered only. [Bottom panel:] the AGN 
duty cycle as a function of the redshift for two galaxy stellar 
masses (logMH,=11.25, 11.75 Mq ). Filled circles and trian- 
gles are previous determination by Brusa et al. 2009 for the 
same masses. The black dashed (dotted) curve is the expec- 
tation of the Menci et al.(2006,2008) SAM for logMH.=11.25 
Mq (logM,=10 Mq ). 



on a stacking analysis of faint X-ray sources in both CDFS 
and CDFN. 

Our results put some first constraints on the faint end of 
the AGN luminosity function (42.75<logL(2-10keV)<44.5) 
at z= 3 — 7, which are already interesting to be com- 
pared to model predictions. To constrain the shape of the 
AGN luminosity function we joined our determinations 
with those obtained at higher luminosities by shallower 
X-ray surveys (Chandra-COSMOS, XMM-COSMOS), and 
optical surveys (SDSS, NOAO DWFS/DLS). We then com- 
pared the broad luminosity range luminosity functions 
with the prediction of the semi-analytic model (SAM) of 
IMenci et al. 20061 IMenci et al. 20081 In this model bary- 



onic processes are associated with the merging histories of 
dark matter (DM) haloes. These haloes contain hot gas at 
virial temperature, a fraction of which can radiatively cool 
down and form a disk with radius and circular veloc- 
ity Vd- During mergers, and also looser galaxy encounters 
(jCavaliere fc Vittorini 2000)) . cold gas can loose its angular 
momentum and can be accreted by the nucleus. A fraction 
of this gas fuels a nuclear starburst while the rest can ac- 
crete on a central SMBH, giving rise to an AGN. The AGN 
timescale is t « r^/wd, i.e. the crossing time for the destabi- 
lized gas. Assuming typical values ^a few kpc and Vd ~ 
100 km/s the AGN timescale is short, a few 10^ yr, com- 
parable with the Salpeter timescale. In this SAM the AGN 
timescale, as well as AGN SMBH masses and Eddington 
ratios, are not free parameters, but are calculated self- 
consistently (the model only assumes that the accretion 
can proceed at most at the Eddington limit). The SAM in- 
cludes a rather detailed tre atment of AGN 'quasar mode' 
feedback (|Menci et al. 2008| . i n terms of a blast wave carry- 
ing outwards the AGN power ( |Lapi et al. 2005[ ). The SAM 
produces prediction for all AGN, both unobscured, mod- 
erately obscured and Compton thick. ILamastra et al. 20101 
compares the SMBH masses and stellar masses of AGN 
host galaxies predicted by this SAM with measurements 
for various galaxy and AGN samples at different redshifts. 
This SAM reproduced reasonably well the z—S-A luminos- 
ity function at logL(2-10keV)> 43.5. At lower luminosities 
it predicts 2-3 times more AGN than observed. It should 
however be noted that some extreme Compton thick AGN, 
those in which the nuclear emission is completely blocked 
by photoelectric absorption and Compton scattering, leav- 
ing only reflection emission in the X-ray band, can well be 
missed even in the deepest X-ray surveys. The agreement is 
sufficiently good at high luminosities (logL(2-10keV)> 45) 
at z=4-5 and z> 5.8 too. In these redshift bins the agree- 
ment between data and model gets worse at lower luminosi- 
ties, where the IMenci et al. 2008) SAM predicts 5-10 times 
more low luminosity AGN than in our determination and 
2-3 times more than the lTreister et al. 2011l ones. 

To investigate the origin of this behaviour and get in- 
sights on high-z AGN physics we also compared our mea- 
surements with predictions derived from more basic models 
for AGN activation through galaxy interactions. This class 
of models consists of two ingredients: a DM halo merger 
rate compatible with cosmological simulations, and an in- 
put AGN light curve (e.g., Wyithe & Loeb 2003, Lapi et 
al. 2006, Shen 2009, Shankar 2009, 2010, and references 
therein). The initial mass of the SMBH at triggering is as- 
sumed to be a fixed small fraction of its mass at the peak of 
activity. SMBH growth is regulated by a condition between 
the peak luminosity and the mass of the host halo at trig- 
gering, consistent with the local relations between SMBHs 
and their host galaxies. The main advantage of approaching 
AGN modeling through this simplified technique is that not 
being part of a specific SAM, they can explore easily and 
quickly a large space of parameters and physical recipes 
to trigger AGN. Here we follow the models presented in 
the preliminary work by Shankar (2010a) with same pa- 
rameters as in Shen (2009). A more comprehensive and de- 
tailed analysis of AGN merger models, with implications 
for SMBH scaling relations and AGN clustering (Shankar 
et al. 2010a, b) is beyond the scope of the present paper and 
will be discussed elsewhere (Shankar et al. 2011, in prepa- 
ration). In this model AGN are activated by mergers of the 
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host DM haloes =M^2/Mft,i > 0.3 for major mergers, 
^ > 0.1 includes minor mergers too). The AGN is trig- 
gered with no dynamical friction time delay between host 
halo and actual galaxy-galaxy merging (a solution close to 
the fly-by hypothesis of lCavaliere fc Vittorini 2000)) . Super- 
Eddington accretion {L/LEdd = 3) is allowed in the initial 
phases of BH growth, and a long, sub-Eddington post-peak 
phase is present in each event. Finally, a minimum halo 
mass Mhmin is assumed (see Shen 2009 for details). We 
plot in the central and lower panel of Fig. [7] four different 
models: 

1. e > 0.3, L/LEdd = 3, M,,™„ = 3 X 10" Afo/Zi (blue, 
solid lines, reference model); 

2. ^ > 0.1, L/LEdd = 1, Mh„^^n = 3 X lO" Mq/Zi (blue, 
dotted lines); 

3. ^ > 0.3, L/LEdd = 3, Mhrmn = 3 X 10" Mq/H, without 
descending phase (blue, thick dashed lines); 

4. e > 0.3, L/LEdd = 1, M,,™„ = 0.1 X 10" Mq/H (blue, 
long-dashed lines); 

The purpose of the comparison of these models with the 
data is to understand whether the merger scenario produces 
expectations in agreement with the data and to highlight 
which of the underlying physical assumptions makes this 
match actually possible. 

The reference model (1) overpredicts by a factor of 2-3 
the AGN luminosity functions at z>4 at all luminosities 
but logL(2-10keV)=43.5-44.5 at z> 5.8. The same model 
without a descending phase (3) reproduces relatively well 
both the z=4-5 and the z> 5.8 luminosity functions. This 
suggests that the underlying AGN light curve must not be 
too prolonged in order not to generate too numerous faint, 
sub-Eddington AGN. 

The inclusion of minor mergers steepens the AGN lumi- 
nosity function at both z—4-5 and z> 5.8. Model (2) with 
^ > 0.1 and L/LEdd = 1 is roughly consistent with the z=4- 
5 luminosity function but underpredicts the bright end of 
the z> 5.8 luminosity function. 

Relaxing the assumption of a minimum halo mass steep- 
ens again the AGN luminosity function at both z=4-5 and 
z> 5.8. The same model with L/Lsdd = 1 (model 4) over- 
predicts the low luminosity end of the z=4-5 luminosity 
function and underpredicts the high luminosity end of the 
z> 5.8 luminosity functions. 

In conclusion, we find that the inclusion of minor merg- 
ers and the exclusion of a minimum halo mass both steepen 
the model luminosity function. About the first, minor 
mergers and galaxy fiy-by helps the formation of stars 
at early times in the progenitors of today massive galax- 
ies, thus helping in explaining their observed colors (e.g. 
IMenci et al. 2004)) . They also help in reproducing quite well 
the X-ray AGN space densities at z ^ 3 (jMenci et al. 2008| . 
About the latter, there could be several causes for a min- 
imum halo mass and/or for inefficiency of BH formation 
and growth in low mass halos. For example, in the SMBH 
formation scenario where seed BHs are formed from the 
monolithic coUaps of gas clouds, the probability for a DM 
halo of hosting a BH seed is an increasing function of the 
halo mass (jVolonteri et al. 2008p . On the other hand, in 
the SMBH scenario where BH seeds are formed from PopHI 
stars, basically ah DM halos with M,, > 10*^"^° Mq at z< 7 
are populated with BHs. Therefore, at least in principle, the 
faint end of the high-z AGN luminosity function might be 



used to discriminate among these two scenarios. However, 
there are at least other three effects that produce a min- 
imum halo mass for BH formation and growth. First, the 
gravitational recoil, giving rise to ejection of one of the BH 
after a merging, can be more important at low masses (see 
e.g. lVolonteri et al. 201ip . Second, a cut-off in the build up 
of galaxies in haloes with Mh < 10^^ Mq could be produced 
by a metallicity effect ([Krumholz &: Dekel 2011|) . Third, an 
intrinsic cut-off at low halo masses can be due to a cut-off 
in the primordial power spectrum of gravitational pertur- 
bation at small scales. Free streaming of warm DM leeds 
to a suppression on scales smaller than the free streaming 
scale and therefore to a cut-off. Warm DM has often been 
invoked to solve problems of the standard ACDM scenarios 
at small scales, such as the large number of galactic satel- 
lites, the cuspiness of galactic cores, and the large number 
of small, compact galaxies predicted by these models (see 
e.g. lPrimack 2009)) . However, the limits on the cut-off scale 
(or, equivalently, on the mass of the warm DM particle) 
derived from high-z Lya forest are quite stringent: up to 
^ 100 kpc no cut-off is seen (jViel et al. 2 008 and refer- 
ences therein). It is difficult to asses quantitatively which 
are the DM halo scales influenced by these effect. To this 
purpose we should include in our SAM and model both 
gravitational recoil and warm DM, and compare their pre- 
diction with AGN and galaxy luminosity functions. This is 
beyond the scope of this work and will be part of a separate 
publication (Menci et al. 2011, in preparation). A minimum 
halo mass for black hole growth implies that accretion is in- 
hibited in halos smaller than this threshold mass (3 x 10^^ 
M0 in our models). BH seeds in these halos do not grow, 
until they are merged to BH in bigger halos. As a conse- 
quence, nuclear accretion occurs only on a fraction of the 
merging history, and the halo population when the above 
threshold is reached late will host SMBH with small masses. 
Qualitatively, they will be somewhat below the Magorrian 
correlation. Again, to obtain quantitative estimates we need 
to study the predictions of models including physical causes 
for the cut-off in the halo mass, which will be addressed in 
a forthcoming pubblication. 

Finally, we have joined our present determination of 
the high-z AGN luminosity function with previous deter- 
mination at lower redshift and with galaxy mass func- 
tion determinations, to empirically estimate the AGN frac- 
tion (or AGN duty cycle), as a function of the galaxy 
stellar mass and redshift. We find that the AGN duty 
cycle increases with both host galaxy stellar mass and 
redshift. The AGN duty cycle is computed by assum- 
ing a given evolution of the Magorrian relationship, 
and by assuming distributions of Eddington ratios at 
different redshifts (see Table 6). These are far to be 
accurately known. Indeed the distributions of Table 6 
taken from[Trakhtcnbrot et al. 2011 > Shemmer et al. 20041 
INetzer fc ^akhtenbrot 2007. Willott et al. 2010bl show an 
evolution with the redshift, while this is not the case 
for the distributions of IKollmeier et al. 20061 However, 
the increase of AGN duty cycle with redshift depends 
marginally on the evolution of the Eddington ratio, refiect- 
ing more basic AGN population properties, evidenced by 
the large evolution of the AGN number density (see e.g. 
IMenci et al. 20061 IMenci et al. 20081 IShankar et al. 2009bl 
and references therein). The correlation of AGN duty cycle 
with both stellar mass and redshift is a rather robust re- 
sult. It holds also by changing the calibration adopted for 



Fiore et et al.: High-z AGN in the CDFS 



17 



the Magorrian relationship by a factor of 2. Our calcula- 
tion confirms what clustering measurements suggest (e.g., 
[Shanka r et al. 201 0a) : the large clustering signal from lumi- 
nous quasars implies that BHs reside in massive haloes that 
by definition are rare and thus require high duty cycles. 

We compare in Fig. ^ the evolution of the AGN duty 
cycle with the prediction of the IMenci et al. 2 008 SAM. 
We see that this model reproduce the observed trend of 
increasing duty cycle from to 5, although the ob- 
served AGN duty cycle at each redshift is somewhat smaller 
than predicted for the same stellar mass. We recall again 
that Compton thick AGN are present in the model, but a 
fraction may be missed in the data. This rough agreement 
suggests that the triggering of at least relatively luminous 
AGN (logL(2-10keV) 2:43 ergs/s), is probably produced by 
galaxy interactions in the galaxy stellar mass and redshift 
range investigated (up to z~ 4 — 5). At higher redshifts the 
model predicts a saturation of the AGN duty cycle of galax- 
ies with masses logM^, > 11.25 Mq (i.e basically all massive 
galaxies at z 5 are expected to host an active nucleus with 
logL(2-10keV)> 43 ergs/s). 



5.1. Future perspectives 

The faintest sources used in the comoving space density 
calculation are from the Chandra-GOODS-ERS sample. 
Indeed, we found that the HST/WFG3 ERS H band images 
are deep enough to trace high-z AGN populations at the ex- 
tremely faint X-ray flux limits reached by the Chandra 4 
Msec exposure, avoiding significant incompleteness, or at 
least with an incompleteness comparable to that reached 
at much brighter fiuxes by Chandra-COSMOS and XMM- 
COSMOS. This is not the case for the HST/ACS optical 
and Spitzer IRAC images from which the GOODS-MUSIC 
catalog is generated. These images are too shallow to find 
the counterparts of all faint X-ray sources at the extremely 
deep CDFS 4 Msec flux limit. We conclude that to fully ex- 
ploit ultradeep Chandra data, i.e. to keep the identification 
rate of faint X-ray sources comparable to that reached at 
brighter fluxes, a deep NIR (R^ 26) coverage is mandatory. 
To avoid large incompleteness in the GOODS area with- 
out IIST/WFC3 coverage we were forced to use only X-ray 
sources brighter than 5 x 10^^*" erg/cm^/s . It must be 
noted that the ERS area lies at relatively large off-axis an- 
gles in the Chandra field, where the sensitivity of Chandra 
is reduced by a factor of at least 2 by both vignetting and 
point spread function (PSF) degradation. As a result, we 
could not exploit in this work the full sensitivity of the 
4 Msec Chandra field. This limitation can be overcome in 
future, when the HST/WFC3 CANDELS survey (Grogin 
et al. 2011; Koekemoer et al. 2011) will be ultimated. The 
CANDELS survey covers the central part of the Chandra 
field and therefore can allow us to fully exploit the Chandra 
sensitivity. Based on the ERS results we expect to at least 
double the sample of faint high-z AGN using the CANDELS 
survey of the CDFS. 

The CANDELS multy-cycle HST treasury program 
covers a total of five fields: the CANDELS deep fields 
are in the GOODS-South and GOODS-North regions, 
and cover a total of 130 arcmin^ to a depth of 11=27.8. 
The CANDELS wide survey covers the AEGIS, COSMOS 
and UDS fields for a total of 670 arcmin^ to a depth of 
H'^ 26.5 (similar or slightly shallower than in the ERS 



field). All these five fields are going to be premiere fields 
for high-z galaxy and AGN studies. If we extrapolate our 
finding on the ERS field to aU CANDELS fields, assuming 
a Chandra coverage of 4 Msec for the CANDELS deep 
fields and of ~ 1 Msec for the CANDELS wide fields, 
we expect 30-100 z=4-5 AGN, 10-30 z=5-6 AGN and 
3-10 z> 6 AGN. This goal can be reasonably achieved 
by Chandra in the next few years (i.e. adding 5 — 6 
Msec of observation to 8 Msec already spent on these 
field). Probably the biggest problem with this sample 
would be the spectroscopic identification of these faint, 
high-z AGN. A fraction can be spectroscopically confirmed 
by HST/WFC3 grism spectroscopy. For the faintest, the 
only viable chances are probably ALMA observations in 
a few years, and, at the end of the decade, JWST and 
ELTs observations. At the end of this program we should 
have a very good knowledge of the faint end of the AGN 
luminosity function up to z=5, but a still rough deter- 
mination of the luminosity function at z=5-7. To make 
further progresses with Chandra, i.e. quantitatively probe 
the first generation of accreting SMBH, which would allow 
putting stringent constraints on SMBH formation mod- 
els pZadau fc Rees 20011 |Lodato fc Natarajan 2006 
Lodato fc Natarajan 2007[ 
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Volonteri fc Begelman 201 0*, 'Begelma n 2010[), and accre- 
tion scenarios (jVolonterf fc Rees 20051 IDotti et al. 20101 
IFanidakis et al. 20111 [King et al. 2008| , would require 
to at least triple the exposure times, i.e. spend on deep 
surveys other 30-40 Msec. While this is not technically un- 
feasible, it is certainly extremely expensive. The Chandra 
limiting problem is that its sensitivity is very good on 
axis, but degrades quickly at off-axis angles higher than 
a few arcmin, making difficult and expensive in terms of 
exposure time to cover with good sensitivity area larger 
than a few hundreds arcmin^. A significant leap forward in 
the field would then be obtained by an instrument capable 
of reaching the Chandra Msec exposure, on axis sensitivity 
(i.e flux limits in the range 1 - 3 X 10-1'^ erg/cmVs ) 
but on a factor of > 10 wider field of view. The extremely 
good Chandra PSF (HPD^ 0.5 arcsec), which is the main 
reason for the good Chandra sensitivity for point sources, 
is beyond reach of today technology of thin X-ray mirrors. 
However, this can be supplied for by increased throughput, 
lower background and wider field of view with nearly 
constant PSF. A straw-man calculation shows that a 
telescope with the total XMM throughput (i.e. area '^half 
square meter) and a PSF ^ 10 times worse than Chandra 
(i.e. HPD~ 5 arcsec) but slowly degrading with the off-axis 
angle, feeding a detector with half degree side and 2-3 time 
lower internal background would be as fast as Chandra but 
with a ~ 20 times better Grasp, i.e. sensitivity times field 
of view. A further improvement could be achieved reaching 
an internal background similar to that of Swift XRT on 
low Earth orbit ( > 10 times smaller than the Chandra 
internal background per unit solid angle). 

The situation at the bright end of the luminosity func- 
tion is fortunately easier. On this side, large areas will be 
surveyed in a few years by the XMM XXL survey and by 
eROSITA. The XMM XXL survey will cover - 25 deg^ 
of sky at a depth of ^ 5 x 10~^^ erg/cm^/s with 10 ks 
effective exposure. eROSITA will survey most of the sky 
down to a 0.5-2 keV fiux limit ~ 10"^'' erg/cm^/s , and 
hundreds deg^ 2-3 times deeper. The identification of high- 
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z AGN in these surveys requires near infrared follow-up. 
This will be possible with both present (VISTA) and future 
(Euclid) instrumentation. In particular, Euclid will provide 
H=24 photometry on most of the high galactic sky and 
H=26 photometry on about 50 deg^. All these surveys will 
help in probing luminous AGN and test the completeness 
of the optical large area surveys such as SDSS. They will 
extremely useful to find rare high-z, high luminosity but 
highly obscured QSOs. 

The present data can only put rough constraints on the 
high-z AGN duty cycle. We stress that the AGN duty cycle 
determinations are presented here for two main purposes: 
first we would like to point out that this measure as a func- 
tion of the redshift is, at least in principle, a powerful tool to 
discriminate among different AGN triggering mechanisms. 
The present data show a fast increase of AGN duty cycle 
in massive galaxy with redshift, naively favoring galaxy en- 
counters as triggering mechanisms. However, the uncertain- 
ties are still too large to draw any strong conclusion. Then, 
we use the AGN duty cycle determination to emphasize 
which are the main areas of improvement to reduce these 
uncertainties. First, we clearly need a better constrain on 
the AGN logL(2-10keV)=43-44.5 density at z> 4. A major 
step toward this direction can be obtained by using further 
combined Chandra and HST surveys in the next few years. 
Second, we need better constraints on the galaxy stellar 
mass functions at z> 3, in particular the low mass end of 
the mass functions. Efforts in this direction are on going, 
using both ground based surveys (with VISTA and other 
near infrared cameras on 8m class telescopes) and from the 
space (HST/WFC3 surveys and the Spitzer SEDS survey). 
Third, we need better constraints on the Eddington ratio 
distributions as a function of the redshift, in particular on 
the distribution at low L/hEdd- Last, we need to test ad- 
ditional physics, such as BH gravitational recoil, and BH 
formation and growth in a warm DM scenario, through the 
comparison of model predictions with observed luminos- 
ity functions and AGN duty cycle. On going and future 
Chandra and HST surveys, in particular those planned in 
the CANDELS fields, should be able to find statistically 
sound samples of z=>5 and logL(2-10keV)> 43 AGN (a 
few tens), and put reliable constraints on the high-z galaxy 
stellar mass function, providing in turn much better con- 
straints on the AGN duty cycle. Fig. IH] also shows the pre- 
diction of the lMenci et al. 20061 iMenci et al.lOOSl SAM for 
low galaxy stellar masses (M* ~ 10^° M©). For such masses 
the duty cycle increase with the redshift, reach a maximum 
and then decreases at z>4. This is due to the combina- 
tion of two effects. First, small galaxies form in relatively 
low density regions, where gas cooling is slower, forming 
smaller gas reservoirs than for large galaxies. Second, at 
the highest redshifts the rate of galaxy interaction is small, 
it increases up to the redshift were groups forms and de- 
crease again at low redshift because high relative velocities 
and low densities. The detection of this distinctive expecta- 
tion would be the unambiguous confirmation of the galaxy 
interaction scenario for AGN triggering. While Chandra is 
sensitive enough to probe ^ 10^ Mq SMBH up to 7, the 
low mass end of the galaxy stellar mass functions at z>4 
is today practically unconstrained. The latter is therefore 
the area where improvements are more needed, but unfor- 
tunately also more difficult with the present generation of 
instruments. A breakthrough in this field would probably 



await for the next generation of infrared telescope such as 
JWST and the ELTs. 



6. Summary 

We have exploited the ultra deep Chandra and HST cov- 
erage of the CDFS to search for X-ray emission at the po- 
sition of high redshift galaxies. To this purpose, we have 
used ephot, a highly optimized tool for X-ray source de- 
tection and photometry (see the Appendix). We found sig- 
nificant emission at the position of 17 galaxies in the area 
covered by the HST/WFC3 ERS survey, and 41 galaxies in 
the full GOODS area for a total of 54 independent z> 3 X- 
ray sources, 29 of which are not present in previous X-ray 
catalogs. We reach a flux limit of ~ 1.7 x lO"'^^ erg/cm^/s 
in the 0.5-2 keV band, which corresponds to luminosities of 
^ lO''^ ergs/s at z=3 and ~ 10^'^ ergs/s at z=7. The present 
Chandra deepest exposure is thus able to probe Seyfert-like 
galaxies up to the epoch of formation of the first galaxies 
and normal star-forming galaxies up to z=3-4. Indeed we 
have evidence that the X-ray emission of one of the faintest 
Chandra-GOODS-MUSIC source, M4417, is likely of stel- 
lar origin. This is the source in our sample with the lowest 
X-ray to H and z bands flux ratio, and the X-ray lumi- 
nosity expected from the SFR rate measured from the UV 
rest emission, oxygen lines and radio 1.4 GHz emission is 
close to the observed one. This can be the highest redshift 
inactive galaxy observed in X-rays. 

We analyzed the X-ray spectra of the Chandra-GOODS- 
ERS sources finding hard spectra in several cases. In 3 cases 
the spectra can be fitted with both a reflection model or 
an absorption model with a column density in excess 10^^ 
cm~^. In an additional case the column density results in 
the range 10^3 - 10^4. The fraction of Compton thick AGN 
at 0.5-2 keV fluxes between 0.3 and 3xl0~^^ erg/cm^/s 
and at z>3 is thus at least 19>t}iQ%- This is higher than pre- 
dicted by popular models for the C XB (Gilli et air^OOT)) . 
Similar conclusions are reached by IGilli et al. 20111 and 
ITreister et al. 20TT1 The optical counterparts of the highly 
obscured and Compton thick AGN do not show any red- 
dening, and we thus conclude that the size of the X-ray 
absorber is likely smaller than the dust sublimation radius 
(^ Ipc, at the luminosity of our AGN). 

We compared the X-ray spectra of the Chandra- 
GOODS-ERS and Chandra-GOODS-MUSIC sources with 
their optical spectra and optical/NIR morphology. Two 
Compton thick or highly obscured AGN show broad ab- 
sorption lines in their UV rest frame spectra. Four of the 
ten highly obscured and Compton thick AGN have a point 
like morphology in their z and H bands images, thus sug- 
gesting that the active nucleus contributes significantly to 
the UV and optical rest-frame emission. Indeed, their rest- 
frame UV (0.16/xm) extinction, estimated by fitting the 
observed SED with galaxy and AGN templates turns out 
smaller than 1-2 magnitudes in all cases. This implies that 
the X-ray absorber must be practically dust free, a con- 
dition satisfied if the X-ray absorbing matter is within or 
close to the dust sublimation radius. This is a fraction of 
pc for our sources, a region more compact than the pc-scale 
obscuring torus. Compact absorbers have been discovered 
in local Seyferts (e.g. IRisaliti et al. 20051) . they might be 
more common at high redshift. 
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We searched for radio emission at the position of 
the Chandra-GOODS-ERS and Chandra-GOODS-MUSIC 
sources in the deep VLA-CDFS maps at 1.4 GHz. We 
find highly significant emission (SNR>5) in 2 Chandra- 
GOODS-ERS sources and 5 Chandra-GOODS-MUSIG 
sources. We find radio emission with lower signal to noise 
(2.4<SNR< 3.2) in additional 8 sources in Table 3. The 
presence of significant radio emission from the faint X-ray 
sources is further confirmed by the analysis of the radio im- 
ages obtained by stacking together the radio images at the 
position of Chandra-GOODS-ERS and Chandra-GOODS- 
MUSIC sources. The average f lux per source of 6-10 /xJy is 
higher than that expected by ILa Franca et al. 20111 based 
on the X-ray fluxes and assuming common nuclear origin. 
This would suggest either a strong ev olution of the AGN ra- 
dio "loudness" with the redshift (the ILa Franca et al. 20111 
determination is mostly based on z< 3 AGN), or that our 
high-z, X-ray detected galaxies, are actively forming stars 
(SFR>a few hundreds Mq/yt). In two cases the SFR eval- 
uated from the radio flux is consistent with that evalu- 
ated from sub-mm observations (M208) or UV SED fitting 
(M4417). 

We evaluated the comoving space density of the z> 3 
faint X-ray sources in three redshift bins: 3-4, 4-5 and 
> 5.8. We found that the wide luminosity range z> 3 lu- 
minosity functions are consistent with a pure luminosity 
evolution model, with L* evolving quite rapidly, and re- 
ducing by a factor 3.3 from z=3 to z=6: L*{z = 3) = 
6.6 X lO^'' ergs/s, L*{z = 6) 2 x 10^4 ergs/s. We 
compared these space densities with the predictions of 
the lMenci et al. 20061 fMenci et al. 2008 SAM and with the 
expectations of more basic models for AGN activation 
through mergers. We find that these models are broadly 
able to reproduce the high-z AGN luminosity functions. A 
better agreement is found by assuming a minimum halo 
mass. We speculate that there are at least four effects or 
scenarios that can produce a minimum halo mass for BH 
formation and growth: 1) if BH are formed from the mono- 
lithic collaps of gas clouds, the probability for a DM halo 
of hosting a BH seed is an increasing function of the halo 
mass; 2) gravitational recoil, giving rise to ejection of one 
of the BH after a merging, which can be more important 
at low masses; 3) metallicity dependent star-formation; 4) 
a cut-off in the primordial power spectrum of gravitational 
perturbation at small scales. To assess the relative impor- 
tance of all these effects we need both better data and more 
complete models for galaxy and AGN formation and evo- 
lution. 

Finally, we empirically evaluated the evolution of the 
AGN duty cycle (fraction of AGN more luminous than 10'*'^ 
ergs/s in the 2-10 keV band, to the total number of galax- 
ies with the same mass) as a function of both host galaxy 
stellar mass and redshift. We compared it with the expec- 
tation of a model using galaxy interaction as AGN trig- 
gering mechanism. The uncertainties on the AGN duty cy- 
cle are large at all redshifts. At z< 4 they are dominated 
by the uncertainty on the Eddington ratio distributions. 
At higher redshift the uncertainty on both AGN luminos- 
ity function and stellar mass functions are extremely large, 
in particular at low AGN luminosities and small galaxy 
stellar masses. Nevertheless, it is clear how the AGN duty 
cycle must change substantially from z=0-0.5 to z=4-5, in- 
creasing by 1.5 orders of magnitude for host galaxy stellar 
masses 1 — 5 x 10^^ Mq. This big rate of change of the 



AGN duty cycle of massive galaxies is roughly consistent 
with that predicted by the models using galaxy interac- 
tions as AGN triggering mechanism. We recall again that 
this applies to relative luminous AGN (logL(2-10keV) J; 43 
ergs/s). Triggering of lower luminosity AGN is of course not 
constrained by our analysis. 

Appendix A: ePhot: energy and space adaptive 
filter for source detection 

The standard techniques for source detections identify 
statistically significant brightness enhancements, deriving 
from both unresolved (point) and resolved (extended) X- 
ray sources, in images accumulated in given energy ranges 
and observation times (such as tools 'celldetect', vptdetect' 
and 'wavdetect' in CIAO, 'emldetect' in SAS). However, 
binning and projection always result in a loss of infor- 
mation. Conversely, present X-ray data, and in particular 
Chandra/ACIS data, are incredibly rich, containing spatial, 
energy and timing information at very good resolutions. 
Even though X-ray data are organized in event files, pre- 
serving the information from the datacube, the scientific 
analysis is usually not performed directly on the original 
event files. The datacube is usually projected onto spatial 
plane before performing source detection, and spectra are 
accumulated at given times or at given sources positions be- 
fore performing spectral fits. To take full advantages of large 
datacubes we should be able to detect sources in a mul- 
tidimensional space. We are developing multidimensional 
source detection techniques using spatial, spectral and tim- 
ing information. We search for clustering of X-ray events in 
a multidimensional space (energy and time in addition to 
the usual spatial coordinates). This allows us to efficiently 
detect X-ray bursty sources and X-ray sources character- 
ized by a peaked spectrum (strong lines, cut-offed spec- 
tra), by reducing the elapsed time and/or the energy range 
(and therefore the background) where the source detection 
is performed. 

We start by accumulating spectra at each galaxy posi- 
tion using different circular source extraction regions from 
3 to 9 pixels on cleaned event files. 

As first step we selected all galaxies with more that 12 
background subtracted counts in the wide 0.3-4 keV band, 
to exclude the faintest objects, and achieve a high reliabil- 
ity (see Sect. A. 2 below). At this stage the background is 
evaluated from a map obtained from the original image af- 
ter the exclusion of all bright sources (sources detected by 
a wavelet algorithm at a conservative threshold). 

For all these galaxies we search for the source extraction 
region radius and contiguous energy band that minimize the 
Poisson probability for background fiuctuation and maxi- 
mize the signal to noise ratio. We adopt a logarithmically 
spaced energy grid with width 0.025 dex from 0.3 to 6.5 
keV. 

Background at the position of each source is evaluated 
by normalizing an average background at the source off-axis 
to the source spectrum accumulated between 7 and 11 keV, 
where the contribution of faint sources is negligible, due to 
the sharp decrease of the mirror effective area. We do not 
use a 'local' background for two main reasons: 1) minimize 
the systematic error due to sharp background variations 
and 2) search for source detections in relatively narrow en- 
ergy bands, where the 'local' background counts may be 
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Table A.l. Chandra background lines 



Energy 


Width 


identification 


keV 


keV 




0.585 


0.05 


Blend of Fe-L and hnes 


1.49 


0.015 


1.4867 Al Ka(+1.5575 Al K^) 


1.78 


0.010 


1.7400 Si Ka (+1.8359 Si K/3) 


2.16 


0.045 


An Map 


7.48 


0.05 


7.47 Ni Ka 


8.26 


0.16 


8.28 Ni KP 


9.71 


0.06 


9.71 Au La 


11.5 


0.1 


11.44 Au L/3 


13.6-14 


0.2 


13.38 Au L7 



CDFS 4MfiCC background s 



small and the associated statistical uncertainty compara- 
tively high. 

A.l. Background evaluation 

The study of faint X-ray sources requires the best possi- 
ble characterization of the background. The typical level 
of the CDFS background in the 0.5-7 keV band is 2-3 
counts/arcsec^ or '^25-40 counts in an area of 2 arcsec ra- 
dius. The spectrum of the Chandra background is com- 
plex, with several broad and narrow components, both in- 
ternal and of cosmic origin. To characterize the Chandra 
background we extracted its spectrum from the following 
regions: 1) circular region with radius 8 arcmin; 2) circu- 
lar region with radius 3 arcmin; 3) annulus with internal 
radius 3 arcmin and external radius 6 arcmin; 4) annulus 
with internal radius 6 arcmin and external radius 8 arcmin. 
Circular region of 5 arcsec radius around all X-ray detected 
sources have been excluded from the analysis in all cases. 
Fig. lA.ll shows the spectrum extracted from the broad 8 
arcmin radius region. The continuum above a few keV is 
dominated by particle induced background (PIB), which 
has a power law shape with positive spectral index above 
4-5 keV. Over-imposed to this continuum we see at least 
8 strong emission lines. The best fit energies of the lines 
are and their identification is given in Table A.l. The ris- 
ing background continuum and the strong emission lines 
present above ^ 7 keV, coupled with the strongly decreas- 
ing mirror effective area above this energy, suggest to limit 
the analysis of faint sources to the band below 7 keV. 

The low energy background spectrum is peaked at 0.5 
keV and decreases rather sharply between ~ 0.6 and ~ 1 
keV. This continuum and the 0.585 line suggest that this 
part of the background spectrum is dominated by the ther- 
mal emission from the local "superbubble" . Below 0.3 the 
response of the instrument decreases sharply and therefore 
we limit the analysis at energies above 0.3 keV. 

We have then fitted the background spectrum between 
0.3 and 7 keV with a model including: 1) three gaussian 
lines to account for the features at 1.48, 1.74 and 2.16 keV; 
2) one broad gaussian line to reproduce the broad bump 
between 1 and 3 keV; 3) two power laws modified at low 
energy by photoelectric absorption to reproduce the broad 
band hard continuum. All these components are not con- 
voluted with the mirror effective area. 4) a thermal compo- 
nent with abundances fixed to Solar value, convoluted with 
the mirror effective area. This model provides a good fit to 
the background spectrum, without systematic residuals. We 
have then fitted the same model to the background spectra 
accumulated in three different regions, adjusting only the 



0.2 0.5 1 2 5 10 

channel energy (keV) 

Fig. A.l. The spectrum of the background in the Chandra 
4 Msec exposure of the CDFS. The spectrum has been ex- 
tracted from a circle of 8 arcmin radius, after the exclusion 
of 5 arcsec radius region around all detected sources (10-20 
arcsec for bright sources and extended sources). 

global normalizations of the model spectra accounting for 
the internal and cosmic background components. The fits 
are remarkably good above 1 keV. Below this energy we 
see a deficit of counts in the 3 arcmin radius spectrum and 
in the 3-6 arcmin radius annulus spectrum, while there is 
an excess in the 6-8 arcmin annulus spectrum. The spectra 
become increasingly soft at large off-axis angles. The rel- 
ative intensity of the thermal component is higher in the 
on-axis spectrum than in the off-axis spectrum, as expected 
if this component is truly cosmic, because of the larger mir- 
ror vignetting at high off-axis angles. The increase of the 
strength of the low energy cut-off at smaller off-axis is then 
probably due to a thicker deposition of contaminants at the 
centre of the ACIS-I 4 chip region. 

The behavior observed for the emission lines is quite 
interesting. While the intensity of the 1.49 keV and 2.16 
keV lines is constant with the off-axis angle, the intensity 
of the 1.78 keV Si line is larger at large off-axis angles. 

The variations of the background spectra accumulated 
in thinner annuli becomes comparable with the statistical 
uncertainty. Guided by this analysis we decided to adopt 
the average background spectra accumulated in the three 
r < 3', 3' < r < 6' and 6' < r < 8' in the following analysis. 

To evaluate the background counts in each source ex- 
traction region we need to normalize the average back- 
ground at three off-axis angles to that at each source posi- 
tion. To this purpose we first extract spectra from circular 
regions of 10 arcsec radii centered at each galaxy positions, 
and then normalize the average background to that in these 
spectra using the 7-11 keV faint-source- free band. The sta- 
tistical uncertainty on the background evaluation is there- 
fore mostly given by the 7-11 keV total number of counts 
at the source position. The choice of the size of the extrac- 
tion region at each galaxy position must then satisfy two 
competing requirements. On one hand the region must be 
large enough to provide enough counts to keep statistical 
fluctuation small. On the other hand, a small region has 
the advantage of avoiding systematic errors due to sharp 
variations of the exposure time (and so of the background) 
typical of mosaics like the ACIS-I CDFS exposures, cos- 
metic defects, crowded regions with several sources within 
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few arcsec. The average CDFS 7-11 kcV background in 
a 10 arcsec radius region is between 600 and 700 counts, 
~ 3 times higher than the 0.5-2 keV average background, 
4 times the 1-2 keV background and 2 times the 1-4 keV 
background. Therefore, a circular region of 10 arcsec radius 
is a good compromise, because it collects a number of 7-11 
keV counts similar or higher than the 0.5-2, 1-2 or 1-4 keV 
counts in 10-20 arcsec annulus (a typical region used to ex- 
tract a 'local' background in crowded areas like the CDFS), 
in a three time smaller area. 

A.2. Completeness, reliability and sky coverage 

To extract statistically quantitative information from a 
given sample of "detected" sources we need to first assess 
its completeness and reliability (number of spurious detec- 
tions). To this purpose we performed a series of extensive 
detection runs on simulated data. 

First, to assess the reliability of our survey we performed 
simulations using background spectra only. We simulated 
about 10^ spectra at the positions of the sources detected in 
the CDFS to use exactly the same exposure time, vignetting 
and PSF, but including only the average background at 
each off-axis position. We run our ePhot on these simula- 
tions and studied the number of spurious detections as a 
function of various parameters: 1) the threshold on back- 
ground subtracted counts in a wide (0.3-4 keV) band; 2) 
the minimum and maximum source extraction regions; 3) 
the energy grid used to search for detections; 4) the band- 
width Emax / Emin'i 5) the Poissou probability threshold for 
a source detection. We choose the combination of parame- 
ters that keep the number of spurious detections ^ 1 every 
5000 spectra, i.e. the number of candidates in our search 
(see Sect. 2): 0.3-4 keV background subtracted counts > 12, 
Emax I Emin > 2, Poisson probability < 2 x 10^. 

To study the completeness of our survey we performed 
detection runs on simulated source plus background spec- 
tra. We used the Mock X-ray catalogs generated by R. Gilli, 
including about 10^ AGN with a distribution of luminosi- 
ties, redshifts and absorbing column densities in agreement 
with the IGilh et al. 20071 CXB model. In particular, the 
Mock catalogs include sources with unabsorbed flux down 
to 10"^® cgs, where the AGN density is - 20000 deg^. 
Source spectra are distributed over six templates, with col- 
umn densities N^f up to 3 x 10^'' cm~^. A pure reflection 
spectrum is also considered to account for heavily Compton 
thick sources. The distributions of spectra is a strong func- 
tion of both luminosity and redshift (see IGilli et al. 20071 
for details). Using these Mock catalogs we simulated a 10 
deg^ Chandra survey with the same characteristics of the 
CDFS survey. We simulated sources at the position of the 
sources detected in the CDFS. We added to each simu- 
lated spectrum a background proportional to the average 
background at each off-axis angle. Finally, we run our on 
each source -I- background simulated spectrum our source 
detection algorithm ePhot. These simulations are used to 
study the completeness of our survey and compute the cor- 
responding sky coverage. Fig. IA.2[ top panel, compares the 
0.5-2 keV observed flux histogram of the simulated sources 
with that of the sources detected by ePhot with the thresh- 
olds defined above. The ratio of the two histograms gives 
the completeness of our survey, or the sky coverage, when 
normalized to the CDFS searched area. Fig. IA.21 bottom 
panel, compares the ePhot sky coverage for the ERS area 
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Fig. A. 2. [Top panel:] the 0.5-2 keV fiux histogram of 
the simulated sources (black solid histogram) and of the 
simulated sources detected by e-phot. [Bottom panel:] the 
Chandra sky coverage in the ERS area obtained using a 
standard detection algorithm (black solid curve) and using 
e-phot (red dashed curve) at the same probability thresh- 
old. 



of the CDFS to that of a traditional detection algorithm at 
the same probability threshold (same fraction of spurious 
detection). 

A. 3. Photometry 

e-phot was also run on the galaxy samples with fixed en- 
ergy bands (0.5-2 keV, thus optimizing only for the size of 
the source extraction region). X-ray fiuxes in the band 0.5- 
2 keV are estimated from e-phot count rates in the 0.5-2 
keV band if the signal to noise ratio in this band is higher 
than 2.5 or from the count rates in the band that optimize 
the detection otherwise. A simple power law spectrum with 
as = 0.4 has been assumed for the count rate flux conver- 
sion. To test our photometry we compared the 0.5-2 keV 
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flux obtained from the event files of the first 2 Msec ob- 
servations to those in Luo08, and the 0.5-2 keV flux ob- 
tained form the full 4 Msec observations to those in Xuell. 
The agreement is remarkable, the median of the two sam- 
ples agree within 3% and the semi-interquartile range of 
their ratio is 0.06. Fluxes of faint detections can be statis- 
tically over-estimated because of the so-called Eddington 
bias dHogg fc TuFileTlQQS, .Wang 2004, ). We estimated the 
bias according to Wang 2004 and found that it is at most 
~ 20% for our faintest detection and it is negligible for 
most of the Chandra-GOODS-ERS and Chandra-GOODS- 
MUSIC sources. Luminosities in Table 3 are calculated from 
fluxes corrected for the Eddington bias. 
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